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Abstract 

We consider asymptotics of the ground state energy of heavy 
atoms and molecules in the self-generatedl magnetic field and derive 
it including Scott, and in some cases even Schwinger and Dirac 
corrections (if magnetic field is not too strong). We also consider 
related topics: an excessive negative charge, ionization energy and 
excessive positive charge when atoms can still bind into molecules. 

This preprint supersedes [Ivr7], [Ivr8], and [Ivr9] as it contains 
more detailed analysis, more results, stronger results and also some 
corrections. 
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Introduction 

We are going to replace Schrodinger operator without magnetic field as in 
Chapter 24 or with a constant magnetic field as in Chapter 25 of [Ivrll] by 
Schrodinger operator 

(0.0.1) H= H Ay = ({D-A)-a) 2 - V(x) 

with unknown magnetic field A but then to add to the ground state energy 
of the atom (or molecule) the energy of magnetic field (see selected term in 



2 



(0.0.2) thus arriving to 



(0.0.2) E(A) = inf Spec(H AV ) + or 1 / |V x A\ 2 dx 



with A/-particle quantum Hamiltonian Y\a,v defined by (25.1.1) of [Ivrll] 
and a parameter a G (0, /t*Z _1 ] with small constant k* > 0. 
Then finally 

(0.0.3) E* = inf EM) 

defines a ground state energy with a self-generated magnetic field 1 ^. 

First of all we are lacking so far a semiclassical local theory and we 
are developing it in Section 1 where we consider one-particle quantum 
Hamiltonian 

(0.0.4) H = H Ay = ({hD-A)-a) 2 - V{x) 

but instead of inf Spec(H,4 y) we consider T^(Ha v) which as we already 
know is what replaces inf Spec(H/\ y) if electrons do not interact (then if 
electrons interact we will need to replace \/ by 1/1/ which includes a potential 
generated by the electron cloud and justify this by estimating an error). We 
define energy of magnetic field as in (0.0.2) but with k replaced by nh d ~ l 
(here d > 2 is arbitrary) and we prove that for d = 2 , 3 in this framework 
self-generated magnetic field is weak and the asymptotics with the remainder 
0(h 2 ~ d ) (or even o(h 1 ~ d ) under standard assumption of the global nature) 
is exactly as for k = (i.e. with 4 = 0). Under standard assumption 
about trajectories we can upgrade this asymptotics to even sharper with 
the remainder estimate o(h~ 1 ) and with term x 2 /? _1 (Schwinger correction 
term) . 

Then in Section 2 we consider operator with potential having Coulomb- 
type singularities and combining results and arguments of Sections 1 and 12.5 
of [Ivrll] prove for d = 3 that 

(0.0.5) lr{H Ay ) = Wey^ + 25(k)/T 2 + 0(k| log + /T 1 ) 

x ) This notion was introduced in series of papers L. Erdos, S. Fournais and J. P. 
Solovej [EFS1, EFS2, EFS3] ; see also L. Erdos and J. P. Solovej [ES3] . 
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provided k < k* (which is a small constant) and there is just one singularity; 
when there are several singularities with a minimal distance a ^> 1 between 
them we prove that 

(0.0.6) Tr{H A _ v ) =\Ney\ 1 + 2S{K)h- 2 + 0(K\\ogK\h-l + h- 1 + K3- 3 h- 2 ). 

As k <C A?5 j log/?| _ 5 under standard assumption about trajectories we can 
upgrade this asymptotics to even sharper with the remainder estimate o(h^ 1 ) 
and with Schwinger correction term. 

Further, in Section 3 we apply these results to provide estimates from 
above and below for the total energy (0.0.3). As a byproduct we also estimate 
D(pvj/ — p TF , pvj/ — p TF ) where W is a ground state for a near-minimizer A. 

This estimate enables us in Section 4 to derive upper estimates for 
excessive negative charge, estimates or asymptotics for ionization energy, 
and in the free nuclei model also for minimal distance between nuclei and 
(in the case of molecule) for excessive positive charge. 



1 Local semiclassical trace asymptotics 

1.1 Toy- mo del 

1.1.1 Statement of the problem 

Let us consider operator (0.0.1) in M d with d = 3 where A, V are real- valued 
functions and V e «5f 2 n J2? 4 , A G Mq. Then this operator is self-adjoint. 
We are interested in Tr~{H A V ) (the sum of all negative eigenvalues of this 
operator). Let 

(1.1.1) E* := inf EM), 

(1.1.2) E(A) := (ir H AV + K^h- 2 J \dA\ 2 dx^j 

with dA = (djAj) a matrix; here and below r is a parameter and constants 
do not depend on it. 
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The estimate from above is delivered by A = and Weyl formula with 
an error 0(/7 _1 ) as V G 2) 

(1.1.3) E* < Weyl-L + 0(h^ 1 ); 

where 

(1.1.4) Weyl(r) = J-h' 3 [{V + tA dx, 

5tt z J 

and 



(1.1.5) Weyl 1= / rd T Weyl(r) = 2 — [ v\dx. 

J-oo ' 157r J 

Also for estimates o(h~ 2 ) we need to include into Wey^ the corresponding 
boundary term. Now our goal is to provide an estimate from below 

(1.1.6) E* > Wey^ - 0(^ 1 ); 

We will use also Weyl(x, r) and Weyl^x) defined the same way albeit without 
integration with respect to x. 

1.1.2 Preliminary analysis 

So, let us estimate E(/\) from below. First we need the following really 
simple 

Proposition 1.1.1. Let V E _Sf i H^ 4 . Then 

(1.1.7) E*>-C/T 3 

and either 

(1.1.8) \: [ \dA\ 2 dx < Ch- 3 



nh 2 



or E(A) > ch' 3 . 



2 - 1 Recall that this means that the second derivatives of V are continuous with the 
continuity modulus 0(| log |x — y|| _1 ), see Section 4.5 of [Ivrll]. If there is a boundary it 
does not pose any problem provided it is in the classically forbidden region. 
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Proof. Using the Magnetic Lieb-Thirring inequality (5) of E. H. Lieb, 
M. Loss, M. and J. P. Solovej [LLS]) 



(1.1.9) J tre 1 (x,x,r)dx > 

- Ch- 3 J Vjdx-Ch 2 J (V 2 J \dA\ 2 dx^(h-* J v A + d X y 

we conclude that for any 5 > 

(1.1.10) E(A) > -Ch' 3 - C5 3 h- 3 + (k- 1 - r 1 )/?- 1 J \dA\ 2 dx 

which implies both statements of the Proposition. □ 
Proposition 1.1.2. Let V + e JSf § n k < ch' 1 and 

(1.1.11) V< -C-\l + \x\) s +C. 

Then there exists a minimizer A. 



Proof. Consider a minimizing sequence Aj. Without any loss of the generality 
one can assume that Aj — > A^ weakly in J^ 1 and in J2? 6 and strongly in 
J2£ c with any p < 6 3 ). Then A^ is a minimizer. 

Really, due to (1.1.10) negative spectra of H Aj V are discrete and the 
number of negative eigenvalues is bounded by A//,. Consider ordered eigen- 
values Aj k of Ha. v. Without any loss of the generality one can assume that 
Xj k have limits A^ < (we go to the subsequence if needed). 

We claim that Aqo./c are also eigenvalues and if Aoo.fc = ... = Aoo,fc+r-i then 
it is eigenvalue of at least multiplicity r. Indeed, let Uj^ be corresponding 
eigenf unctions, orthonormal in J2? 2 . Then in virtue of Aj being bounded in 
Jzf 5 and V G =5f 4 we can estimate 

\\Du Jik \\ < KWuj.kWl" 7 ■ \\uj,k\\ a < KWDujA 1 -* ■ \\uj, k \\° 

with a > which implies HDt/y^H < K. Also assumption (1.1.11) implies 
that ||(1 + IxD^Uy /tH are bounded and therefore without any loss of the 
generality one can assume that Uj^ converge strongly. 

3 ) Otherwise we select a converging subsequence. 
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Then 

(1.1.12) lim Tr{H A . v ) > Tr{H Aoo y), 

(1.1.13) liminf J \dAj\ 2 dx > J |<9A»| 2 dx 

and therefore E(/4oo) < E*. Then A^ is a minimizer and there are equalities 
in (1.1.12)— (1.1.13) and, in particular, there no negative eigenvalues of H AooV 
other than A^. □ 

Remark 1.1.3. We do not know if the minimizer is unique. Also we do not 
impose here any restrictions on r, K (which may depend on h) in (1.1.11) or 
k > 0. From now on until further notice let A be a minimizer. 

Proposition 1.1.4. In the framework of Proposition 1.1.2 let A be a mini- 
mizer. Then 

(1-1-14) ^A^-(x) = ; := 

-tr^2(vjV k (hD k - A k ) x + G k Oj{hD k - A k )^je{x,y, r)\ y=x 

k 

where A = (Ai, A 2 , A 3 ), a = (<Ti, cr 2 , o" 3 ) and e(x, y, r) is the Schwartz kernel 
of the spectral projector Q(—H) of H = H AV and tr is a matrix trace. 

Proof. Consider variation 8A of A and variation of Tr"(H) = Tr(/-/~) where 
H = HQ( — H) is a negative part of H. Note that the spectral projector of 
His 



;i.l.l5) 9(r - H) = -L T Res R (r - H)' 1 



and therefore 



5Tr(9(r-H)) = -^ f Res R Tr((r - H)-\bH){r - H)' 1 ) = 

- [ T Res R Tr((8H)(r- H)- 2 ) = -d T ^-. f Res K Tr((6H)(r - H)' 1 ) 
^' J—oo Ziyi j— oo 



- d T Tr ({5H)9{t - H)). 
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Plugging it into 

/o po 
rd T Tr(e{r - H)) = - Tr(0(r — H) dr) 
-oo J — oo 

and integrating with respect to r we arrive after simple calculations to 

(1.1.17) STr-(H) = Tr((6H)9(r- H)) = -J <$>{x)8A{x) dx 
where ^>{x) is the right-hand expression of (1.1.14). Therefore 

(1.1.18) 8E(A) = [ (-<D(x) - -^-AA(x))8A{x) dx 

J nn z 

which implies (1.1.14). □ 
Proposition 1.1.5. If for k = k* 

(1.1.19) E* > Weyli - CM 
mt/j M > C/r 1 then for k < k*(1 - e ) 

(1.1.20) \ [ \dA\ 2 dx< CiM. 

«n J 

Proof. Proof is obvious based also on the upper estimate E* < Wey^ + 
Or 1 . □ 

1.1.3 Estimates 

Proposition 1.1.6. Let estimate (1.1.20) be fulfilled and let 

(1.1.21) <,=nMh<c. 
Then as r < c 

(i) Operator norm in Jz? 2 of (hD) k d(r — H) does not exceed C for k = 0, 1, 2; 

(ii) Operator norm in J2? 2 of (hD) k ((hD — A) ■ cr)0(r — H) does not exceed 
C fork = 0,1. 



Proof, (i) Let u = 9(r - H)f. Then ||u|| < ||f || and since 

(1.1.22) \\A\\x6 < C\\dA\\ < C(kM)H 
we conclude that 

\\hDu\\ < \\(hD - A)u\\ + \\Au\\ < \\{hD - A)u\\ + C\\Ay 6 ■ \\u\\s* < 
\\{hD - A)u\\ + C{ K M)h\\u\\^ 2 ■ \\u\\f 6 < 
\\{hD-A)u\\ + C{KMh)^\\u\\ 1/2 ■ \\hDu\\ 1/2 < 

\\{hD - A)u\\ + ^\\hDu\\ + CnMh\\u\\; 

therefore due to (1.1.21) 

(1.1.23) \\hDu\\ < 2\\{hD - A)u\\ + CKMh\\u\\. 
On the other hand, for B = V x A and r < c 

\\(hD - A)u\\ 2 < C\\u\\ 2 + (h\B\u, u) < C\\u\\ 2 + h\\B\\ ■ \\u\\%, < 
C\\u\\ 2 + C(KM)h 2 \\u\\ ■ \\uye < C\\u\\ 2 + C{KM)h\\u\\ ■ \\hDu\\ < 

C(l + nMh 2 + JMh 2 )\\u\\ 2 + ^\\{hD - A)u\\ 2 

and due to (1.1.23) we conclude that 

(1.1.24) \\{hD - A)u\\ < C\\u\\ and \\hDu\\ < C(l + KMh)\\u\\ 

provided nMh 1+s < c for sufficiently small S > 0. Therefore under assump- 
tion (1.1.21) for k = 0, 1 statement (i) is proven. 

Further, since (hD) 2 = (hD - A) 2 + A(hD - A) + AhD - h[D, A] we in 
the same way as before (and using (1.1.24)) conclude that 

\\{hD) 2 u\\ < C\\u\\ 2 + ^\\hD(hD - A)u\\ + ^||/? 2 D 2 ty|| 

and therefore 

\\h 2 D 2 u\\ < C\\u\\ 2 + C\\AhDu\\ 

and repeating the same arguments we get ||/7 2 D 2 u|| < C||u||; so for k = 2 
statement (i) is also proven. 

(ii) Statement (ii) is proven in the same way. □ 
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Corollary 1.1.7. Let (1.1.20) and (1.1.21) be fulfilled. Then as r < c 

(1.1.25) e(x,x,r) < C/T 3 

and 

(1.1.26) |((/7D-/\)-(j)e(x,y,r)| x=y | < C/T 3 . 

Proof. Due to proposition 1.1.6 operator norms from Jzf 2 to ^ of both 
0(r — H) and ((/7D — /4) • a)0(r — H) do not exceed C and the same is true 
for an adjoint operator which imply both claims. □ 

Corollary 1.1.8. Let (1.1.20) and (1.1.21) be fulfilled and A be a minimizer. 
Then 

(1.1.27) \\dA\\vi- S < CkH- 1 
and 

(1.1.28) \\dA\y~ < C' 5 h-t- 5 

where is the scale of Holder spaces and 5 > is arbitrarily small. 

Proof. Really, due to (1.1.14) minimizer A satisfies || A/4||j?oo < Cnbr 1 . Also 
we know that \dA\ < C{KMh 2 )\ < Ch\ due to (1.1.21). Then (1.1.27) holds 
due to the standard properties of the elliptic equations 4 - 1 . 

Therefore if at some point y we have |9/4(y)| > //, it is true in its 
e(/i/?/t~ 1 ) 1 ~ <5 -vicinity (provided \i < /t/?" 1 ) and then 

and we conclude that 

fi 2 ^- 1 ) 3 ^ < Cnh 2 M « fi 5 ~ 3S < CK 4 - 3S h- 1+3S M 

and one can see easily that (1.1.28) holds due to (1.1.21) and assumption 
h- 1 < M< h- 3 . 

On the other hand, if /i > nh' 1 then we need to take e-vicinity and then 
fi 2 < CnMti 2 < Chi where we used (1.1.21) again. Therefore (1.1.28) has 
been proven. □ 

4 ) Actually we can slightly improve this statement. 
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Remark 1.1.9. (i) It is not clear if it is possible to generalize this theory to 
arbitrary d > 2 with magnetic field energy given by 

(1.1.29) ^ J (\dA\ 2 - \V ■ A\ 2 ) dx 

Surely one should use generalized Pauli matrices cry in the definition of the 
operator: as d = 2 one can prove that E(/\) is bounded from below and 
minimizer exists; as d = 4 one can prove that E(A) is bounded from below 
as k < e h; especially problematic is the case d > 5; 

(ii) Therefore while arguments of Subsection 1.2 below remain valid for 
d > 4, so far they remain conditional (if a minimizer exists and satisfies 
some crude estimates). 



1.2 Microlocal analysis unleashed 
1.2.1 Sharp estimates 

Now we can unleash the full power of microlocal analysis but we need to 
extend it to our framework. It follows by induction from (1.1.27)— (1.1.28) 
and the arguments we used to derive these estimates that 

(1.2.1) \\dA\\ V n- S < CnKh- 1 -" Vn e Z+ 

so A is "smooth" in e = h scale while for rough microlocal analysis as in 
Section 2.3 of [Ivrll] one needs at least e = Ch\ \ogh\. We consider in this 
section arbitrary d > 2; see however remark 1.1.9. 

Proposition 1.2.1. For a commutator of a pseudo-differential operator 
with a smooth symbol and ^ e+1 -function A(x) a usual commutator formula 
holds modulo 0(/? 6,+1 |||9/4||| 6) ) for any non-integer 9 > where 



[1.2.2) 



supl^^x)! #eZ H 

a:\a\=e * 



sup\x-y\ w - ■\d a f{x)-d a f(y)\ 6 £ Z+. 

*-.\ a \=m x ^ y 

Proof. Easy proof is left to the reader. □ 
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Proposition 1.2.2. Assume that 

(1.2.3) \\dV\\ nB ( x ,i)) < Q 

and 

(1-2.4) fi := ||9y4||^ (B(Xjl)) < C . 

Let U(x,y, t) be the Schwartz kernel of e lh ltA . Then for T x 1 

(i) Estimate 

(1.2.5) \F t ^ h -i TX T(t){hD x ) a (hD y fU\ < Ch 5 

holds for all a : \a\ < 2, (3 : \(3\ < 2, s as \x — y\ > C Q T and r < Cq; 

(ii) Estimate 

(1.2.6) \F t ^ TX T{t){hD x y{hD y y Vl {hD x ) V2 {hD y )U\ < 

Ch^Ch e {lA\\ e+1 + \\V\\ e+1 ) 

holds for for all for all a: \a\ < 2, (5 : \/3\ < 2, s and all (pi,<p2 G ^o°°; such 
that dist(supp (pi, supp (pi) > CqT, and t < c ; 

(Hi) If also in in B(x, 1) 

(1.2.7) eo<|^|<c 
then for a small constant T = e estimate 

(1.2.8) \F t ^ T XT(t)(hD x r(hD y fU\<Ch s +Ch e (\lA\\ g+1 + l\Vl\ e+1 ) 

holds for for all for all a : \a\ < 2, : < 2, s as \x — y\ < €qT and 
\t\ < e. 

Proof. Let u = e lth lH f with arbitrary f. 

(i) Statement (i) is easily proven by the same arguments as in the proof of 
Theorem 2.1.2 of [Ivrll]: we consider just usual function <ft(x) and operators 
of multiplication like x(0( x )) so there are no "bad" commutators due to 
non-smoothness of A or V. 
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(ii) Statement (ii) is also proven by the same arguments; however in this case 
= 0(x, £) so we need to involve "bad" commutators but their contributions 
are bounded by 

in the right-hand expression while the left-hand expression is e /7 1 1 u 1 1 2 
where Q, Qi, and Q' are operators with symbols x(0( x >O)' Xi(0( x .O)i an d 
Xi(0( x >O ~~ v) respectively, r\ > is an arbitrarily small constant (so the 
latter symbol has a bit larger support than the former one), 5 > is a small 
exponent, Xi(0 = (~ x'(0)*> an d f.e. x(0 = e ~' f ' 1 as t < 0, x(0 = as 
t > 0. 

Therefore we conclude that 

\\Qu\\ < Ch d (lAj e+1 + \iV\l +1 )\\u\\ + Ch s \\Q'u\\ 

and similarly we can estimate || Q'u\\ with || Q"u\\ in the right-hand expression 
etc and thus we conclude that 

\\Qu\\ < Ch 6 (jA\i 0+1 + \jVi\ e+1 \\u\\ + Ch s \\u\\ 

which is what we need. 

(iii) Statement (iii) is easily proven by the same arguments as in the proof of 
Theorem 2.1.2 of [Ivrll]: we consider just usual function 0(x) and operators 
of multiplication like x(0( x )) so there are no "bad" commutators due to 
non-smoothness of A or V. However we need to consider a contribution of u 
which is not confined to the small vicinity of (y , if) and we need Statement (ii) 
for this so the last term in the right-hand expression of (1.2.8) is inherited. 

We leave easy details to the reader. □ 

Remark 1.2.3. (i) Statement (i) means the finite propagation speed with 
respect to x; 

(ii) Statement (ii) means the finite propagation speed with respect to £ 
and the last term in the right-hand expression of (1.2.6) is due to the 
non-smoothness of A and V; 
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(iii) Statement (iii) means that under assumption (1.2.7) there actually is a 
propagation with respect to x; 

(iv) So far we have not assumed that V is very smooth function; we actually 
do not need it at all: it is sufficient to assume that dV is very smooth in 
microscale e = h 1-6 ; one can actually invoke more delicate arguments of the 
proof of Theorem 2.3.1 of [Ivrll] and deal with microscale e = Ch\ log h\. 

Therefore in the framework of Proposition 1.2.1 (iii) estimate 

(1.2.9) \F t ^ TXT {t){{hD x y{hD y yU{x,y,t))\ x=y \ < 

C^-d+sj-s + CT 2 h- d+s {\\A\l 9+1 + \\V\l e+1 ) 

holds for all a : \a\ < 2, (3 : \/3\ < 2, s as T = e and \r\ < e where as usual 

X6«g B ([-l.-|]U[|.l]) ) XT(t) = x(t/n 

Let us consider T G (Ch, e); then we apply the standard rescaling 

t i — v tT^ 1 , x i — v xT^ 1 , h i — v hT^ 1 and assumptions (1.2.3), (1.2.4) are 

replaced by weaker assumptions 

(1-2.3)' T\\dV\\ nB{Xtl)) < C 

and 

(1.2.4)' T\\dA\\ nB{Xil)) < C . 

Further, ||/4|| e+1 and ||| acquire factor T e+1 . 

Furthermore, as U(x,y, t) is a density with respect to y we need to add 
factor T~ d to the right-hand expression and due to F t _j./,-i T we need to add 
another factor T and after these substitution and multiplications we arrive 
to 

Proposition 1.2.4. Let h < T < e and assumptions (1.2.3)', (1.2.4)' and 
(1.2.7) be fulfilled. Then estimate (1.2.9) holds. 

Next we apply our standard arguments: 
Proposition 1.2.5. In the framework of proposition 1.2.4 

(1.2.10) \F t ^ h - W [x T (t){(hD x r(hD y fU(x,y,t))]\ x=y \ < 

Ch'- d +CT 2 h- d+s (\lA\l d+1 + \\V\l e+1 ) 
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provided x £ ^o°°([ — 1- 1]) an( ^ 

(1.2.11) | [((/,D X - A(x)) ■ a) a ((hD y - A(y)) ■ afe(x, y, t)} \ x=y - 

Weyl ai/S (x)| < 

^(1+ iiMii v(fl(Xil „ + \\av\\ vmxil)) ) + ^-^i^di^ii^+i^ii^,) 1 

u>/iere 

(1.2.12) Weyl^x) := const h~ d [ ((£ - /A(x)) • a) a+P d£ 

J{H(x,Z)<t} 

is the corresponding Weyl expression and 

(1.2.13) H(x, = ((£ - A(x)) -a) 2 - V(x); 

in particular \Ney\ a /3 (x) = as \a\ + \(3\ = 1. 

Proof. Obviously summation of (1.2.9) over CqH < \t\ < T and a trivial 
estimate by Ch 1 ~ d of the contribution of the interval \t\ < CqH implies 
(1.2.10). 

Then the standard Tauberian arguments and (1.2.10) imply that the 
left-hand expression of (1.2.11) does not exceed 

CT-W- d + CTh- d +\lA\\ 9+1 + \\V\\ e+1 ). 

Optimizing with respect to T < e such that (1.2.3)', (1.2.4)' hold we pick 
up T = T* with 

(1.2.14) T* = 

emin(l, (||^||^ ( b ( x,i)) + I|5\/||^(b(x,i))) _1 , T^Q^IIm+l ■ 

Meanwhile the Tauberian formula and (1.2.9) imply that the contribution 
of an interval {t : |t| x 7"} with h < T < T* to the Tauberian expression 
does not exceed the right-hand expression of (1.2.9) divided by T, i.e. 

Ch i-d+s T -s-i + CTh- d+e \\dA\\ v e; 

summation over 7"* := h 1 ^ 5 < T < T* results in the right-hand expression 
of (1.2.11). 

So, we need to calculate only the contribution of {t : \ t\ < T*} but one 
can see easily that modulo indicated error it coincides with \Ney\ a/3 . □ 
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Remark 1.2.6. As d > 3 one can skip assumption (1.2.7). 



Indeed, we can apply the standard rescaling technique: x i— > x£ , 
h h = h£- 3 2, A Af\, V h-> V^ 1 with I = max(e| ihv~\\ 
v = (1 + see Section 5.1 of [Ivrll]. 

1.2.2 Application 

Let us apply developed technique to estimate a minimizer. 
Proposition 1.2.7. Let k < c and let A be a minimizer. Let 

(1.2.15) fi := \\dA\\ v < Ch~ 1+S . 

As d = 2 let assumption (1.2.7) be also fulfilled. Then as 6 £ (1, 2) estimate 

(1.2.16) \\dA\y e -i + h^WdAye < 

C K (l + \\V\\^+h^-^V\\l e+1 ) + C\\dA\\' 

holds with 

(1.2.17) \\9A\\':= sup 11^11^2(6(3,,!)). 

y 

Proof. Consider expression for AA. According to equation (1.1.14) and 
Proposition 1.2.5 \AA\ + \hdAA\ does not exceed the right-hand expression 
of (1.2.11) multiplied by Ck/?^ 1 i.e. 

(1.2.18) \\AA\y + \\hdAA\y < 

Ck(i + \dA\ v + \avy + h^ -v\\dA\\l e + h^ e -v\\dv\\l^ 

where we replaced and by larger ||<9/4||<^e and respec- 

tively. 

Then the regularity theory for elliptic equations implies that 

(1.2.19) For any 9' G (1,2) /7 e ' -1 ||cM||<^0' does not exceed this expression 
(1.2.18) plus C\\dA\\'. 
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Note that ||9/4||^ does not exceed e||<9/4||^e + C e '||<9/4||' with arbitrarily small 
constant e > and therefore 

(1.2.20) h^WdAye+e-'WdAy 

does not exceed expression (1.2.18) plus C e '||d/4||' where we used (1.2.19) for 
6' = 6. 

Comparing (1.2.20) and (1.2.18) we conclude that for k < c and suffi- 
ciently small constant e > we can eliminate in the derived inequality both 
contributions of dA to (1.2.18) thus we arrive to (1.2.16). □ 

Having this strong estimate to A allows us to prove 

Theorem 1.2.8. Let k < c, (1.2.15) be fulfilled, and let d = 3. Assume 
that 

(1.2.21) \\Vye + i<c 
with 9 e (1,2). Then 

(1.2.22) E* = Weyli + 0{h 2 - d ) 

and a minimizer A satisfies 

(1.2.23) \\dA\\ < C«M 
and 

(1.2.24) + tf-^WdAWv* < Cith* + Ck. 

Proof, (a) In virtue of (1.2.10) the Tauberian error when calculating Tr(H^ v ) 
does not exceed the right-hand expression of (1.2.10) multiplied by CT~ 2 
i.e. 

(1.2.25) Ch^T' 2 + Ch- d +°(m\e + i + IMIU)- 

Assumption (1.2.21) allows us to simplify this expression and take T x 
(1 + /i) -1 ; applying estimate (1.2.16) we conclude that the Tauberian error 
does not exceed 

(1.2.26) C(l + /2) 2 h 2 ~ d + C(k + \\dA\\')h 2 - d . 
We claim that 
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(1.2.27) Weyl error 5 ) when calculating Tr(H A v ) also does not exceed (1.2.26). 



Then 



(1.2.28) E(A) > 

Wey^ - C(l + fi) 2 h 2 - d - C(« + ||«94||> 2 " d + K- l h l - d \\dA\\ 2 > 



because /x < C||9/4||' + 1 due to (1.2.16) and assumption (1.2.21). This 
implies an estimate of E* from below and combining with the estimate 
E* < E*(0) = Weyli + Ch 2 - d from above we arrive to (1.2.22) and (1.2.23) 
and then (1.2.24) due to (1.2.16) and assumption (1.2.21). 

(b) To prove (1.2.27) let us plug A e instead of A into ei(x,x,0). Then 
in virtue of the rough microlocal analysis contribution to Weyl error of 
{t : 7"* < \t\ < e} with 7"* = h 1 ^ 5 would be negligible and contribution of 
{t : \ t\ < 7"*} would be Wey^ + 0{h 2 - d ). 

(c) Now let us calculate an error which we made plugging A e instead of 
A into ei(x,x,0). Obviously it does not exceed Ch~ d \\A — A £ \\<jg and since 



A - A £ \\ v < Ce e+1 \\dA\\^e this error does not exceed Ch e+1 - d - A5 \\dA\\^e 



which is marginally worse than what we are looking for. 

However it is good enough to recover a weaker version of (1.2.22) and 
(1.2.23) with an extra factor h~ Sl in their right-hand expressions. Then 
(1.2.16) implies a bit weaker version of (1.2.24) and in particular that its 
left-hand expression does not exceed C. 

Knowing this let us consider the two term approximation. With the above 
knowledge one can prove easily that the error in two term approximation 
does not exceed Ch 3 - d - 5 ' with 5' = 1005. 

Then the second term in the Tauberian expression is 



where subscript ( £ ) means that we plugged A £ instead of A and superscript 
T means that we consider Tauberian expression with T = T* = e. But 




(1.2.29) 




5 ) I.e. error when we replace Tauberian expression by Weyl expression. 



18 



then the contribution of {t : T* < \ t\ < T*} is also negligible and modulo 
Ch e+2 ~ d ~* s \\dA\\%>e we get a Weyl expression. However 

(1.2.30) (H AV - H AeM ) = -2(£ - A £ ) ■ (A - A £ ) + \A- A £ \ 2 

and the first term kills Weyl expression as an integrand is odd with respect 
to (£ — A £ ) while the second as one can see easily makes it smaller than 
Ch 3 ~ d ~ s '. Therefore (1.2.27) has been proven. □ 

Remark 1.2.9. (i) For d = 2 we cannot drop assumption (1.2.7) at the stage 
we did it for d > 3. However results of the next section allow us to cure this 
problem using partition-and-rescaling technique. 

(ii) Actually as V G c € 2,1 we have an estimate 

(1.2.31) \dA(x) - dA(y)\ < Ck\x - y|(| log |x - y\\ + 1) + C/i. 
Combining with (1.2.23) we conclude that 

(1.2.32) \\dA\\ v < C^ d+1 ^ d+2 ^\ log/ ? | d /( d+2 )/) 1 /( d+2 ) 

(hi) If (1.2.21) holds for (hd) m V with m e Z+ then (1.2.23) and (1.2.23) 
also hold for (hd) m A instead of A; further, if (hd) m V e then (1.2.31) 
and (1.2.32) also hold for (hd) m A instead of A. 

1.2.3 Classical dynamics and sharper estimates 

Now we want to improve remainder estimate 0(h 2 ~ d ) to o(h 2 ~ d ). Sure, 
we need to impose condition to the classical dynamical system and as 
\dA\ = 0(h s ) with 5 > due to (1.2.32) it should be dynamical system 
associated with the Hamiltonian flow generated by Hoy: 

(1.2.33) The set of periodic points of the dynamical system associated with 
Hamiltonian flow generated by H v has measure on the energy level 0. 

Recall that on : H Q y(x,£,) = r} a natural density d\x T = dxd£ : 

dH\ H=T is defined. 

The problem is we do not have a quantum propagation theory for H A V as 
A is not a "rough" function. However it is rather regular function, almost ^ 2 , 
and (A-A e ) is rather small: \A-A £ \ < r] := Ch 2 ~ 35 and \d(A-A £ )\ < Ch l ~ 35 
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and therefore we can apply a method of successive approximations with the 
unperturbed operator Ha e ,v as long as r]T/h < \f i.e. as T < h 1 ^ 46 . Here 
we however have no use for such large T and consider T = 0(h~ s ). 
Consider 

(1.2.34) F t ^ h -i T x T (t)U(x,y,t), 

and consider terms of successive approximations. Then if we forget about 
microhyperbolicity arguments the first term will be 0(h~ d T), the second 
0{h- l - d r]T 2 ) = Oih 1 -*- 5 ') and the error 0{h- 2 - d r] 2 T 3 ) = 0{h 2 - d - 5 "). 

Therefore as our goal is 0(h 1 ~ d ) we need to consider the first two terms 
only. The first term is the same expression (1.2.34) with U replaced by U^ £ y 

Consider the second term, it corresponds to U'^(x,y, t) which is the 
Schwartz kernel of operator 

(1.2.35) U' (e) := ih- 1 f e^-^ h ~ lHA - v (H AV - e''^ 1 "^ dt' 

Jo 

and then 

(1.2.36) Tr(U' (£) V>) = //T 1 Tr((H AV - H A , U ) e''"^ ^(t)) 
with 

JO 

is /)-pseudo-differential operator with a rough symbol and V'HO ~ t. 

Really, one can prove easily studying first the Hamiltonian flow equation 
and then the transport equations that ip t := e lh tH ^yipe~ lh tHA e* is a 
/7-pseudo-differential operator with a rough symbol and its corresponding 
norm is bounded. 

Note that 

r 1 F t ^- lT e' rl ^^ 1 (t) = (20 J {F t ^ T ,e ih - ltH ^)f{h-\r - r')) dr' 

with f = F t ^ T f t , f t = xtW 1 ^) and therefore (1.2.35)-(1.2.36) imply that 

(1.2.37) |F t ^,- lTXT (t)TrU' (£) V>| < C V T 2 h~ d 

where in comparison with the trivial estimate we gained factor h. 

We can plug here 7"' e ( T* , T) instead of T and taking summation by 
V from 7"* = e to T we conclude that (1.2.37) also holds for Xr{t) replaced 
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by (xr(t) — Xt*(0) (provided x = 1 on (— |, |)) and since 77 T 2 < /? 1+<5 as 
7" < the right-hand expression (1.2.37) does not exceed Ch 1 ~ d+S . 
On the other hand, our traditional methods imply that as d > 3 

(1.2.38) |F t ^ / ,- lTXT (t)Tr(e /t/ '" 1 ^.^)l < Ch l ' d 7>(n T , c ) + C T ^h l - d+5 

where Y\j is the set of points on energy level 0, periodic with periods not 
exceeding T, \^t,c is its (- vicinity, ( > is arbitrarily small. 

Here again we can plug any T' e ( 7"* , T) instead of T and after sum- 
mation with respect to V we conclude that (1.2.38) also holds with Xr(t) 
replaced by (xr(0 - Xr»(t))- 

Combining with estimate for (e' th Ha ^ m — e' th HA e^ we conclude that 

|F t ^- lT (xr(0 - xr.(t)) Trfe*'*- 1 ^^)! < Ch^ 7>(n T , c ) + C^/, 1 ^ 
and since 



we conclude that 



(1.2.39) \F t ^ h - lTXT (t)Jr(e it,h - lH -^)\ < 

Ch x - d + Ch 1 - d Tyi{n TX ) + Cr^/) 1 -^ 

and then the Tauberian error does not exceed the right-hand expression of 
(1.2.39) multiplied by ChJ- 2 and it is less than CT^h 2 ~ d . 

Consider now the Tauberian expression and again apply two-term ap- 
proximation for e lh th>A - v considering e lh tHA ? v as an unperturbed operator; 
then the error will be less than Ch 2 ~ d+S . 

Consider the second term after taking trace; it is 0(h 2 ~ d ~ 4S ), so it is just 
slightly too large. Further, if ip = I one can calculate it easily and observe 
that it is 0{h 2 - d+5 ) provided V e < ^ 2 ' 1 . 

Finally, the first term is what we get for e lh tH ^y and in virtue of 
rough microlocal analysis contribution of {t : 7~* < |t| < T} does not 
exceed Ch 2 ^ d \i{V\ T ,^) + C T ,c,h 2 ~ d+& and contribution of {t : |t| < T*} is 
Weyli + 0{h 2 ~ d+s )'. 

Then we arrive to 
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Theorem 1.2.10. Let k < c, (1.2.15) and (1.2.21) be fulfilled, and let 
d = 3. Furthermore, let condition (1.2.33) be fulfilled 6 **. Then 



(1.2.40) 



E* = WeylJ + o{h 2 - d ) 



where 



(1.2.41) 



Weyl* = Weyl x + x/? ; 



I 



VlAVdx 



d 



calculated in the standard way for Ho v and a minimizer A satisfies similarly 
improved versions of (1.2.23) and (1.2.24). 

Remark 1.2.11. (i) Under stronger assumptions to the Hamiltonian flow 
one can recover better estimates like 0(h 2 ~ d \ \og h\~ 2 ) or even 0(h 2+s ~ d ) 
(like in Subsubsection 4.4.4.3 "Sharper remainder estimates") of [Ivrll]. 

(ii) We leave to the reader to calculate the numerical constants x* here and 
in (1.3.5) below, x = K\ — |x 2 . 

(iii) However, even if ip ^ I we can observe that it is sufficient to consider 
only principal terms and then the second term in approximations is also 
0(h 2 ~ d+s ) provided V G c € 2 ' y as long as principal symbol of ip(x) is even 
with respect to £, in particular, if ip = ip(x). 

1.3 Local theory 

1.3.1 Localization and estimate from above 

The results of the previous Subsection have two shortcomings: first, they 
impose the excessive initial requirement (1.1.21) to k as a priory M < ch~ 3 ; 
second, they are not local. However curing the second shortcoming we make 
the way to addressing the first one as well using the partition and rescaling 
technique. 

We can localize Tr (H) = Tr(H 1 ) which is the first term in E(A) ei- 
ther in our traditional way as Tr(H~?p 2 ) or in the way favored by some 
mathematical physicists 7 -* : namely, we take Tr~(ipHip) where in both cases 
ip G < ^ oo (6(0, !)), <ip <1 and some other conditions will be imposed to 



it later. 



6 ) I.e. no(noo) = 0. 

7 > See f. e. L. Erdos, S. Fournais. and J. P. Solovej [EFS1]. 
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Note that 



(1.3.1) Tr-(^H^) > Tr{ipH-iP) = J ei (x, x, 0)^ 2 (x) cbc. 

Really, decompose operator H = HQ(-H) + H(l - Q{-H)) where 6(r - H) 
is a spectral projector of H and therefore in the operator sense H > H~ : = 
HQ( — H) and ^> /-/■?/> > ifjH~if) and therefore all negative eigenvalues of 
are greater than or equal to eigenvalues of the negative operator ipH~i/j and 
then 

r° 

(1.3.2) Tr^HiP) > Tr(ip H~ ip) = Tr / tc! t Q{t - H)ip 2 

which is exactly the right-hand expression of (1.3.1). 
Remark 1.3.1. Each approach has its own advantages. 

(i) In particular, no need to localize A (see (ii)) and the fact that proposi- 
tion 1.1.5 obviously remains true are advantages of Tr~(-?/;H^)-localization. 

(ii) Further, as Tr~(ipHip) does not depend on A outside of 6(0, |) we may 
assume that A = outside of 6(0, 1). Really, we can always subtract a 
constant from A without affecting traces and also cut-off A outside of 6(0, 1) 
in a way such that A' = A in 6(0, |) and \\dA'\\ < c||<9/4||e(o,i); the price is 
to multiply n by c _1 -as long as principal parts of asymptotics coincide. 

(iii) On the other hand, additivity rather than sub-additivity (1.4.1) and 
the trivial estimate from the above are advantages of Tr(-?/> /-/"■?/> )-localization. 

We will use both methods and here we provide an upper estimate 
for larger expression Tr~ (?/>/-/?/>) and a lower estimate for lesser expression 
Tr(ipH~if)). Let us estimate from the above: 

Proposition 1.3.2. Assume that V e c € 2 ' 1 , d > 2. Let £(x) be a scaling 
function 8 ^ 1 and tp be a function such that \d a ip\ < cip£' a ^ for all a : \a\ < 2 
and < cr( 1+5 ) with a>l and 5 > 9 \ 

8 ) I.e. l> and \di\ < \. 

9 ) Such compactly supported functions obviously exist. 
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Then, as A = 0, 

(1.3.3) Tr(^H^) = J Weyl 1 (x)^ 2 (x) dx + 0(/> 2 - d ) 

and under assumption (1.2.33) 

(1.3.4) Tr(^H^) = J Weyl*(x)^ 2 (x) cbc + o(/? 2 ~ d ) 

with 

(1.3.5) Weyl*(x) = Weyl^x) + ^h^V^AV + x 2 h~ 1 Vf l \VV\ 2 

calculated in the standard way for H o v . 

Proof. Let us consider H = ipHip as a Hamiltonian and let e(x,y, r) be the 
Schwartz kernel of its spectral projector. Then 

(1.3.6) Tr(ipHip) = J e 1 {x,x,0)dx = J2 J h{x,x,Q)i)j dx 



where ip 2 form a partition of unity in W 1 and we need to calculate the right 
hand expression. The problem is that H is not a usual Schrodinger operator 
because of degenerating factor ip on each side. 

Consider first an e£-admissible partition of unity in 6(0, 1). Let us 
consider 7-scale in such element where 7 = e£ a and we will use 1 scale 
in £. Then after rescaling x 1— > X7 _1 the semiclassical parameter rescales 
h h-> h new = /?7 _1 and the contribution of each 7-subelement to a semiclassical 
remainder does not exceed Cip 2 (h/ r y) 2 ~ d with ijj < 7 1+<5 having the same 
magnitude over element as 7 > 2h. Then contribution of ^-element to a 
semiclassical error does not exceed Cip 2 (h/ / y) 2 ~ d x£ r y 2 ~ d x Ch 2 ~ d i/j 2, y~ 2 £ d < 

Ch 2-d £ d+25 

Note that expression (1.3.6) only increases if we sum only with respect 
to elements where £ a > h. Therefore we arrive to estimate 

Tr(^H^) < J \Ney\ 1 {x)^j 2 (x)dx+ Ch 2 - d 

where integration is taken over domain {x : £{x) > h 1 ^}. Note that we can 
extend this integral to M. d : really, it will add negative term with absolute 
value not exceeding Ch~ d x h 2+s as ip < h 1+s there and it is absorbed by 
the remainder estimate. □ 
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Corollary 1.3.3. In the framework of proposition 1.3.2 

(1.3.7) E; := inf E^{A) < J Weyl 1 (x)^ 2 (x) dx + Ch 

and under assumption (1.2.33) 

(1.3.8) E; < J Weyl 1 (x)^ 2 (x) dx + Ch 



2-d 



,2-d 



with 



(1.3.9) E^A) := Tr"(^/#) + J \dA\ 2 dx. 
Really, we just pick A = 0. 

1.3.2 Estimate from below 

Now let us estimate redefined E^(A), 

(1.3.10) E4A) := J e 1 ( x ,x,0)^ 2 (x)c/x+-^ ZT J \0A\ 2 dx. 

from below. However we need an equation for an optimizer and it would be 
easier for us to deal with even lesser expression involving r-regularization. 
Let us rewrite the first term in the right-hand expression in the form 

(p(T/L)rd T e(x,x,T)+ / (1 - <^(t/L))t d T e(x, x, r) > 

■oo J — oo 

/ (<f(r/L)(r - L) d T e(x, x, r) + (1 - <p{t/L))t d T e(x, x, r)) 

J —CO 

where (p G ^ 7 oo ([— 1, 1]) equals 1 in [—§,§] and let us estimate from below 

(1.3.11) E;(4):= j(j L y{T/L)(r-L)d T e{x,x,T){x)+ 

(l-^r/L)){T-L)d T e(x,x,T)y\x)dx +1 J^- d j \0A\ 2 dx 
Let us generalize Proposition 1.1.4: 
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Proposition 1.3.4. Let A be a minimizer ofE'^(A). Then 

(1-3-12) J-AAj{x) = *j:= 

tr {^ajO k {hD k - A k ) x + o k Oj{hD k - A k ) y ) x 

/ SK \<p{t/L){t - L) Res K (r - H)^ 2 (r - 

J — oo 

(1 - v(t/L))t(t - L) Res K (r - H)~^ 2 (r - H)" 1 ] (x, y) t/r ) 

J y=x/ 

where we use a temporary notation SK[B](x,y) /or the Schwartz kernel of 
operator B. 

Proof. Follows immediately from the proof of Proposition 1.1.4. □ 

Proposition 1.3.5. Let d = 3 and assumptions (1.1.20) and (1.1.21) be 

fulfilled. Then as r < c 

(i) Operator norm in J?? 2 of (hD) k (r — H) -1 does not exceed C\ lmr| _1 for 
/c = 0,l,2; 

(ii) Operator norm in Jzf 2 of (hD) 2 {fhD — A) ■ cr)(r — H)" 1 does not exceed 
Cllmr]- 1 for k = 0,1, 2. 

Proof. Proof follows the same scheme as the proof of Proposition 1.1.6. □ 

Proposition 1.3.6. Let d = 3 and assumptions (1.1.20) and (1.1.21) be 
fulfilled. Then |<D(x)| < Or 3 . 

Proof. Let us estimate 

(1.3.13) | J t<p(t/ L) Res R SK T(r - H)-^ 2 {t - H)' 1 {x,y)dr\ 

where L < c and </? G ^^([—l, 1]) and also a similar expression with a 
factor (r — /.) instead of r; here either T = I , or T = (/jD^ — /U) x or 
T = {hD k - A k ) y . 

Proposition 1.3.5 implies that the Schwartz kernel of the integrand does 
not exceed Ch~ 3 \ lmr| -2 and therefore expression (1.3.13) does not exceed 
CL 2 x h~ 3 L- 2 = Or 3 . 
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Then what comes out in <t> from the term with the factor (f)(r/h) does 
not exceed Ch~ 3 . 

Then representing (l — 0(r//?)) as a sum of ip{r / L) with L = 2"h with 
n = 0, ... , [I log h\\ + c we estimate the output of each term by Ch~ 3 and 
thus the whole sum by Ch~ 3 \ \ogh\. 

To get rid off the logarithmic factor we use equality 

(1.3.14) (r - H)-V(r - H)' 1 = -d{r- H) V + (r - H)- 2 [h, ^](r - H)" 1 ; 

if we plug only the second part we recover a factor /)//_ where /) comes from 
the commutator and 1/L from the increased singularity; an extra operator 
factor in the commutator is not essential. Then summation over partition 
results in Ch~ 3 . 

Plugging only the first part we do not use the above decomposition but 
an equality Res R (r — H)^ 1 dr = d T d(r — H). □ 

Corollary 1.3.7. Let d = 3, (1.1.20) and (1.1.21) be fulfilled and A be a 
minimizer. Then (1.1.27) and (1.1.28) hold. 

Proof. Proof follows the proof of corollary 1.1.8. □ 

Now we can recover both proposition 1.2.7 and main theorems 1.2.8 
and 1.2.10: 

Theorem 1.3.8. Let d = 3 and assumptions (1.1.20) and k < c be fulfilled. 
Then 

(i) The following estimate holds: 

(1.3.15) E; - J Weyl 1 (x)V' 2 (x) dx = 0{h 2 - d ) 

and and a minimizer A satisfies (1.2.23) and (1.2.24); 

(ii) Furthermore, let assumption (1.2.33) be fulfilled (i.e. (^(rioo) = 0). 
Then 

(1.3.16) E; - J Weyl*(x)^ 2 (x) dx = o(h 2 - d ) 

and a minimizer A satisfies similarly improved versions of (1.2.23) and 
(1.2.24). 
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1.4 Rescaling 

We consider only d = 3 here. 
1.4.1 Case k < 1 

We already have an upper estimate: see Corollary 1.3.3. Let us prove a 
lower estimate 10 ^. Consider an error 

(1.4.1) ( J Weyl^x)^ 2 dx - E^(/\)) + . 
Obviously Tr~ is sub-additive 

(1.4.2) Tr-(^^^y) > ^Tr(^W^) 

and therefore so is E^(/4) under assumption that ^ G ^ 2 (6(x / -, |^-)) where 
multiplicity of covering by B(xj,£j)) does not exceed Co and we are allowed 
to replace k by C\K in the right-hand expression 11 ). 

Then we need to consider each partition element and use a lower estimate 
for it. While considering partition we use so called ISM identity: as 

(1.4.3) Z>y =1 

j 

we have 

(1.4.4) H = J2(i>jH4>j+ 1 -[[H,4>jUj}) = 

j 

j k 

where the second equality is due to the fact that [[/-/, ipj\, ipj is an ordinary 
function. 

10 > But only for E^(A) denned by (1.3.9). 

") Really, ||cM|| 2 > c£. ll^jll 2 with Aj(x) = (A(x)-Aj(xjMj with Vj € S?(S(x;, 
equal 1 in 6(x,-, f •<!/))• 
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In virtue of Proposition 1.1.5, from the very beginning we need to consider 

(1.4.5) M = /^/T 1 -" 

with a = 2, f3 = and k < c. But we need to satisfy precondition (1.1.21) 
which is then 

(1.4.6) n p+1 h~ a < c. 

Therefore, if condition (1.4.6) is fulfilled with a = we conclude that the 
final error is indeed 0(/) _1 ) or even o(/j _1 ) (under assumption (1.2.33)) 
without any precondition. 

Let precondition (1.4.6) fail. Let us use 7-admissible partition of unity 
ipj with ipj satisfying after rescaling assumptions of Proposition 1.3.2. 

Note that rescaling x (—>■ X7 -1 results in h 1— >■ /? new = It-/' 1 and after 
rescaling in the new coordinates ||<9/4|| 2 acquires factor r ) d ~ 2 and thus factor 
K^ 1 h~ 2 becomes k' 1 h~ 2r y d ~ 2 = ^newVew with k h-» /t new = kj. 

Then after rescaling precondition (1.4.6) is satisfied provided before 
rescaling K ^h~ a ^ a+ I 3+1 < c. Let us pick up 7 = «-(/ 3 + 1 )/(«+0+ 1 )/i a /(«+0+ 1 ). 
Obviously if before rescaling condition (1.4.6) fails, then h <^ 7 < 1. 

But then expression (1.4.1) with ip replaced by ipj does not exceed Ch~^ w = 
C(/77 -1 ) -1 and the total expression (1.4.1) does not exceed C(/77~ 1 )~ 1 7^ 3 = 
C/rV 2 = C/^'/r 1 -"' with 

(3' = 2{f3 + l)/(a + (3 + 1), a' = 2a/(a + (3 + 1); 

So, actually we can pick up M with a, (3 replaced by a', f3' and we have 
a precondition (1.4.6) with these new a',/3' and we do not need an old 
precondition. Repeating the rescaling procedure again we derive a proper 
estimate with again weaker precondition etc. 

One can see easily that a' + (3' + 1 = 3 and therefore on each step 
a + (3 + 1 = 3 and we have recurrent relation for a': a' = |a; and 
therefore we have sequence for a which decays and becomes arbitrarily small. 
Therefore precondition (1.4.6) has been reduced to k < h 5 and estimate 
M = 0(/) _1 ) has been established. After this under assumption (1.2.33) we 
can prove even sharper asymptotics. 

To weaken assumption k < h 5 to k < c we can use rescaling x 1— > x^ 1 
with 7 = h s . We arrive to the error estimate 0(h~ 1 ~ s ) and therefore 
optimizer satisfies ||V x A\\ < h^~ s (where 5 is increased if necessary but 
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remains arbitrarily small). Then instead of ||A/4||_^oo = 0(1) we arrive to 
||A/\||^oo = 0(h' 5 ) and to ||9 2 /\||^oo = 0(h~ 5 ); then \\dA\\^ = 0(h^ 5 ); 
it is more than sufficient to unleash microlocal analysis without any need 
to appeal to Proposition 1.1.6 which is the only place where we needed 
assumption (1.1.21). 
Thus we arrive to 

Theorem 1.4.1. Let d = 3, V G ff 2,1 , k < c and let ip satisfy assumption 
of proposition 1.3.2. Then 

(i) Asymptotics (1.3.15) holds; 

(ii) Further, if assumption (1.2.33) is fulfilled then asymptotics (1.3.15) 
holds; 

(Hi) If (1.3.15) or (1.3.16) holds for E^(A) (we need only an estimate from 
below) then \\dA\\ = 0((nhY) or \\dA\\ = o((/t/?)2) respectively. 

1.4.2 Case 1 < k < h 1 

We can consider even the case 1 < k < h" 1 . The simple rescaling-and- 
partition arguments with 7 = k' 1 lead to the following 

(1.4.7) As 1 < k < Z?" 1 remainder estimate 0{n 2 h~ l ) holds and for a 
minimizer A satisfies ||<9/4|| 2 < Cn 3 h. 

However we would like to improve it and, in particular to prove that if k 
is moderately large then the remainder estimate is still O^h" 1 ) and even 
o(/? _1 ) under non-periodicity assumption. 

Theorem 1.4.2. Let d = 3, V G , and let ip satisfy assumptions of 
proposition 1.3.2. Then 

(i) As 

(1.4.8) k < K* h := e/T?| \ogh\-* 
asymptotics (1.3.15) holds; 

(ii) Furthermore as k = o(n* h ) and assumption (1.2.33) is fulfilled then 
asymptotics (1.3.16) holds; 
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(Hi) As K* h < k < ch 1 the following estimate holds: 

(1.4.9) |E; - J Weyl 1 (x)V' 2 (x) dx\ < C/T 3 (/t/))i| log/t/?| 2 

Proof, (i) From (1.2.16) we conclude as k > c that 

h l - e \dA[ge < Ck{k + JI). 

Then using arguments of subsection 1.2.2 one can prove easily that for 
k < \f~\ 

\F t ^H-wXT(t)(hD x nhD x y (U(x,y,t)-U {£) (x,y,t)-U[ £) (x,y,t))\ < Ch 1 '" 

where we use the same 2-term approximation, T = Let us take then 

x = y, multiply by s~ d ip{e~ 1 {y — z)) and integrate over y. Using rough 
microlocal analysis one can prove easily that from both U^(x,y, t) and 
U( £ ){ x < y< f ) we S e t 0{h~ 2 ) and in the end of the day we arrive to the estimate 
\AA e \ < CkJI which implies 

(1.4.10) \d 2 A e \ < CkJi\ log h\ + C/i 

where obviously one can skip the last term. Here we used property of 
the Laplace equation. For our purpose it is much better than |<9 2 /4 e | < 
Ck 2 \ log h\ + C/i which one could derive easily. 

Again using arguments of subsection 1.2.2 one can prove easily that 

(1.4.11) | Tr(ip H A V ip) - Tr(^H- < Qi 2 ^ 
and therefore 

(1.4.12) | Jr{ip H A V ip) - J Weyl 1 (x)^ 2 (x) dx\ < Cfl 2 h 2 - d 

and finally for an optimizer 

(1.4.13) \\dA\\ 2 < CKjl 2 h. 

Here fi and fx were calculated for A, but it does not really matter as due 
to |<9M| < CK 2 h~ s we conclude that \dA - dA £ \ < CK 2 h~ 5 e < C due to 
restriction to k. 
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Then, as d = 3 
(1.4.14) Af 2 (/VM log ^|)) 3 < Kjl 2 h 

and if fi > 1 we have fi = [i and (1.4.14) becomes /t~ 3 | log /7|~ 3 < Cnh which 
impossible under (1.4.8). 

So, n < 1 and (1.4.14) implies (1.3.15) and (1.4.13), (1.4.14) imply that 
for an optimizer ||<9/4|| < C(/t/?)2 and /i < C/t 4 /)| log/)^. So (i) is proven. 

(ii) Proof of (ii) follows then in virtue of arguments of subsection 1.2.2. 

(iii) If n* h < k < /r 1 we apply partition-and-rescaling. So, h i-)- h' = /?7 _1 
and k i— > k,' = k,'-/ and to get into (1.4.8) we need 7 = e/t~3/j~3| log(/t/))| _1 
leading to the remainder estimate Ch~ lr y~ 2 which proves (ii). □ 

Problem 1.4.3. Repeat arguments of Subsubsections 1.1.2.2 and 1.1.3.3 
and of this Subsection as d 7^ 3. When they hold? 



2 Global trace asymptotics in the case of 
Coulomb-like singularities 

2.1 Problem 

We consider the same operator (0.0.4) as before in M? but now we assume 
that V has Coulomb-like singularities. Namely let y m G M 3 (m = 1, ... , M , 
where M is fixed) be singularities ("nuclei"). We assume that 

(2.1.1) V= J2 Z ™ + W ( x ) 

l<m<M l X ~y m l 

where 

(2.1.2) z m >0, Z! + ...+z M xl, 

and 

(2.1.3) \D a W\<C a z m (|x-y m | + l)" 1 |x-y m r H 



l<m<M 

\/a : \a\ < 2 
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but at the first stages we will use some weaker assumptions. Later we assume 
that V(x) decays at infinity sufficiently fast. Let us define E*) and E(A) by 
(1.1.2)-(1.1.1). Finally, let £(x) min 1 < m < M £ m (x) with £ m (x) := ||x-y m |. 
In this and next Sections we assume that 

(2.1.4) re £ (0, re*] where < re* is a small constant. 

As re = we set A = and consider E := Jr~(H A V ); then our results will 
not be new. 

2.2 Estimates of the minimizer 

Let us consider a Hamiltonian with potential V and let A be a minimizing 
expression (1.1.2) magnetic field. We say that A is a minimizer and in the 
framework of our problems we will prove it existence. 

2.2.1 Preliminary analysis 

First we start from the roughest possible estimate: 

Proposition 2.2.1. Let V satisfy (2.1.1)-(2.1.3) and let re < re*. Then the 
near-minimizer A satisfies 

(2.2.1) | J (tr e A i(x, x, 0) - Weyl^x)) dx\ < Ch~ 2 
and 

(2.2.2) \\dA\\<CK*. 

Proof. Definitely (2.2.1)-(2.2.2) follow from the results of L. Erdos, S. Four- 
nais, and J. P. Solovej [EFS3] but we give an independent easier proof based 
on our Subsection 1.1.1. 

(i) First, let us pick up A = and consider Tr(^9(— H)^) with cut-offs 
i/; e (x) = ip((x-y m )/£) where ip G ^{B(0, 1)) and equals 1 in 6(0, \). Here 
and below 9(r — Ha,v) is a spectral projector of H. 
Then 

(2.2.3) | Tr(^H~ u (0)^) | < C7T 2 as £ = L ■= h 2 . 
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On the other hand, contribution of B(x, £) with £(x) > £* to the Weyl 
error does not exceed C(, 2 h~ l = C( 3 £h- 1 where H = h/((£) in the rescaling; 
so after summation over £ > £* we also get 0(h~ 2 ) provided ( 2 < C£~ x . 
Therefore we arrive to the following rather easy inequality: 

(2.2.4) | J (tr e ,i(x, x, 0) - Weyl^x)) dx\ < Ch' 2 . 

This is what rescaling method gives us without careful study of the singu- 
larity. 

(ii) On the other hand, consider A ^ 0. Let us prove first that 

(2.2.5) Tr~ (faHipt) > Ch' 2 - Ch' 2 J \0A\ 2 dx as £ = £*. 

Rescaling x h-> (x — y m )/£ and r h-> rjl and therefore h h-> x 1 

and /\ t— >■ A£^ (because singularity is Coulomb- like) , we arrive to the same 
problem with the same k (in contrast to Subsection 1.4 where k h- > k£ 
because of different scale in r and h) and with £ = h = 1. 

However this estimate follows from L. Erdos, J. P. Solovej [ES3] (we 
reproduce Lemma 2.1 of this paper in Appendix 5. A. 

(iii) Consider now ip£ as in (i) with £ > £*. Then according to theorem 1.4.1 
rescaled 



(2.2.6) Tr (if> e H AiV i/> e ) > -CCHh- 1 - Ch' 2 [ 

JB 



\dA\ 2 dx. 

B(x,2£/3) 



Really, rescaling of the first part is a standard one and in the second part we 
should have in the front of the integral a coefficient K,- 1 h- 2 C 2 x £~ 2 £(h/££)~ 2 
where factor ( 2 comes from the scaling of the spectral parameter, factor 
(- 2 comes from the scaling of the magnitude of A, factor £ = £ 3 x £~ 2 
comes from the scaling of dx and d respectively, and K := h/((£) is a 
semiclassical parameter after rescaling. Therefore this expression acquires a 
factor C 2 ^ < C. 
Then 

(2.2.7) J (tre A i(x,x,0)-Weyl 1 (x))c6c> -Ch~ 2 - Ch 2 J \8A\ 2 dx 
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and adding the magnetic field energy k 1 h 2 \\dA\\ 2 we find out that the 
left-hand expression of (2.2.1) is greater than the same expression with 
A = plus {k- 1 - C)h- 2 \\dA\\ 2 minus C\r 2 which implies (2.2.1) and (2.2.2) 
as A is supposed to be a near-minimizer. □ 

Proposition 2.2.2. Let V satisfy (2.1.1)— (2.1.3). Then there exists a min- 
imizer A. 

Proof. After Proposition 2.2.1 has been proven we just repeat arguments 
of the proof of Proposition 1.1.2. If V e =5f 2 no change would be required 
but for V f. 5£ 2 one needs to consider modifications as in Remark 2.2.3(i) 
below. □ 

Remark 2.2.3. We are a bit ambivalent about convergence of J Weyl-^x) dx 
at infinity, as for Coulomb potential it diverges. In this case however we can 

5 

either assume in addition that V e J2? 5 5 or tackle it as in Proposition 2.4.4 
below. 

2.2.2 Estimates to a minimizer. I 

Let us repeat arguments of Subsection 1.1.3.3. However our task now is 
much more complicated: while we know a priory that ||<9/4|| 2 < Ck we will 
not be able to improve it significantly (or at all as k x 1). 

Recall equation (1.1.14) for an minimizer A. After rescaling x h-> x/£, 
t i->- t/C 2 , h i->- h= h/((£), A h-> A^i this equation becomes 

(2.2.8) AA = 

-2K( 2 m 2 Y,t r (vjOk{hD k -C 1 A k ) x + v k v j {hD k -C 1 A k ) y )e{x,y,T)\ y=x ^ 

k 

and since so far ( 2 £ = 1 we arrive to 

(2.2.9) AA = 

- 2Kh 2 J2 tr (°j°k{frD k - C"M fc ) x + o k Oj(hD k - C 1 A k ) y )e(x, y, r)| y=x ) . 

k 
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(i) Plugging for u = ipd( — H)f with cut-off function if) and repeating argu- 
ments of Subsubsection 1.1.3.3 we conclude that in the rescaled coordinates 

(2.2.10) \\(HD X ■ a)u\\ < \\((HD X - A) ■ a)u\\ + C\\A\\#i ■ \\u\\^ 

< \\{{hD x -A) ■ a)u\\ + Ch-^\\A\\^ ■ \\u\\* ■ \\hD x ii\\* 

< \\{{hD x - A) ■ a)u\\ + ^D x ii|| + C(/H||/\||^ 6 ) 2 ||iv|| 

where ||/4||_^>6 calculated in the rescaled coordinates is equal to || Ax-igll^orig 
(where subscripts "orig" means that the norm is calculated in the original 
coordinates and A) which does not exceed Ck? due to (2.2.2) 12 ) and therefore 
(since \\(hD x ■ a)u\\ = ||ftD x fv||) 

(2.2.11) \\hD x u\\ < C{1 + KfT l )\\f\\. 

Continuing arguments of Subsubsection 1.1.3.3 we conclude that in the 
rescaled coordinates 

(2.2.12) \\{hD x ) k u \\ < c(i + Kh-^wq, 

(2.2.13) \\(hD x ) k {{hD x - A) ■ a)u\\ < C(l + ^fllll- 
for k = 0, 1, 2 and therefore 

(2.2.14) ||A/\||^oc ( e (x ,i)) < CkK~\1 + k/T 1 ) 3 , 
Here we estimate different norms of A locally. Then either 

(2.2.15) ||^||^ (e(x , f)) + ^||9M||^ (B(Xi!)) < CKh-\l + nh- 1 ) 3 
or 

(2.2.16) ||9/\||^=o (s(Xil _ e)) + ^||9 2 /\|| 

< C||(9/4|| = C||(9/4 or ig||ori g < Ck^ 
In the latter case (2.2.16) we have in the original coordinates 

(2.2.17) 11^11^=0(8^) < C«iri 

12 ) As usual we assume that the average of A over B(x, 1) is 0. 
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and we are rather happy because then the effective intensity of the magnetic 
field in B(x,£) is C"^ll^||^=o (B{Xil _ e)) < CkI 

(ii) The former case (2.2.15) is much more complicated because our estimate 
is really poor as k x 1 and we are going to act only in this assumption. 
Assume that 

(2.2.18) ||9/4||j5foo(e( x ,i-e)) < H 

with fj, > fr a . Selecting u = ipd(—H)f with 7-admissible ip we conclude 
that \\(A ■ cf)u\\ < ||/4||jfoo||iy|| < C/ryH u|| (assuming without any loss of the 
generality that A = at some point of suppVO and that 

\\{hD) k u\\ < C(l + %- 1 + /i7)\ 
\\{hD) k {{hD - A) ■ &)u\\ < C(l + + ^) k+1 

and therefore 

\V x (hD x - A) • tr)e(., ., 0)| < Cr 3 (l + + ni) 1 

and then 

||A^||^oo( S ( Xi i_ eH < Ch~\i + h^ 1 + ni)*. 

Optimizing with respect to 7 = yT^\rh we conclude that either 

II^IU^x,!-,)) < Ch-^ii + htf 

or (1.1.21) holds. In the former case using the second of estimates 

(2.2.19) ||4||*» {fl(Xil _ e)) < C\\d 2 A\\^ {Bixl _ £)) \\dA\\l, 

(2.2.20) ||^|U~(B(x,i- e) ) < C^Ali^^JdAfs 
we conclude that (2.2.18) with 



and one can see easily that starting from /j, = h 4 as given by (2.2.15) we 

3 j~ 

can arrive after number of iterations to fi = h~$ and therefore 
(2.2.21) ||9M||^oc (s(Xil _ e)) < Ch-^ 1+2 V~ S k = 0,l, 2. 
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(iii) This estimate is good enough to launch microlocal arguments. Assuming 
(2.2.18) with /i < \f^ x we estimate as in Section 1 

|r x ((ftD x -/\).(j)e(.,.,0)| < CfifT* 

and then 

\\d 2 A\\ B{xA _ e) < Cfih~ s 

and therefore 

||^|| e(x ,i_ e) < Cf,h- 5 

resulting in fx := fx^h~ 5 and after a number of iterations we get fx = h~ 5 and 
therefore iterating this procedure one more time and taking into account 
factor k we arrive to 

(2.2.22) Either (2.2.16) holds or 

(2.2.23) II^IUoc^,!^)) < CKh~ s . 

However to prove that the effective magnetic field 0(1) we need to modify 
these arguments, and we do it in the next Subsubsection. 

2.2.3 Estimates to a minimizer. II 

In this step we repeat arguments of Subsubsection 1.2.1.1 but now we have 
a problem: we cannot use [i = ||cM||oo as we have domains X r = {x : 
£(x) > r} rather than the whole space. So we get the following analogue of 
(1.2.18) where A is still rescaled and the norms are calculated in the rescaled 
coordinates: 

(2.2.24) llA^ll^aj + HWAdAW^s^ < 

C«(l + |cMk( e(x ,i) + ^^ll^ll^e^)) 

which implies 

(2.2.25) \\dA\\ v{B{Xil)) + H^WdAW^s^ < 

eh^- l \- l \\dA\\, g o {B{xA)) + CK\\dA\\ nB(x , 1)} + C||9/\||^2 (s(Xil)) 
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and the last term in the right-hand expression does not exceed Cxk. 

Let u(r) = sup x .^ x ) >r f(x) where f(x) is the left-hand expression of 
(2.2.11) calculated for given x in the rescaled coordinates. Then (2.2.25) 
implies that for k e (0, k*) (where k* > is a small constant) 

v{r) < \v{\r) + C^ 

which in turn implies that 

v{r) <^v{2- n r) + 2C^, n > 1, 

and therefore 

z/(r) < 4Cks + 4 sup f(x) < Qk^ 

C h 2 <£(x)<2C h 2 

due to the rough estimate (because h x 1 as £(x) x h 2 ). Then returning 
to the original (not rescaled) coordinates and to the original (not rescaled) 
potential A we arrive to estimates (2.2.26) and (2.2.27) below: 

Proposition 2.2.4. Let k < k* , ( = ci~^ . Let A be a minimizer. Then 
for £(x) > 4, = h 2 estimate (2.2.17) holds and also 

(2.2.26) \d 2 A(x) -d 2 A(y)\ < CkU~\ \x - y\ £ 0/2 £; e/2 < 6 < 1, 
and 

(2.2.27) \dA(x) - dA(y)\ < CkSt** \x - y|(l + | log |x - y\\). 
Remark 2.2.5. (i) So far we used only assumption that 

(2.2.28) \d a V\ < CC 2 r H Va : \a\ < 2 
with ( = but even this was excessive; 

(ii) In this framework however we cannot prove better estimates as (2.2.17) 
always remains a valid alternative even if ( <^ £'2^ 

(iii) Originally we need an assumption (1.2.7) \ V\ > e , but for d = 3 one 
can easily get rid off it by the standard rescaling technique. 

Consider now zone {x : £{x) < 

Proposition 2.2.6. Let k < k* , ( < c£~^. Let A be a minimizer. Then 
\dA\ < CkH~ 3 as £{x)<£* = h 2 . 
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Proof. Proof is standard, based on the rescaling (then h = 1) and equation 
(1.1.14) for a minimizer A. We leave easy details to the reader. □ 

Let us slightly improve estimate to a minimizer A. We already know that 
\dA(x)\ < C f3 with p = £—2 and using the standard rescaling technique we 
conclude that 

(2.2.29) \AA\ < Ck( 2 (3 + CnC^t 1 

5 

which does not exceed Ck£~2 which implies 
Proposition 2.2.7. In our framework 

(i) As £{x) > h 2 

(2.2.30) \A\ < CkTK \dA\ < &d~* 
and 

(2.2.31) \dA{x) - 8A{y)\ < Qk^-^x - y\ e as \x - y\ < ^£{x) 
for any 9 G (0, 1); 

(ii) As £(x) < £* = h 2 these estimates hold with £(x) replaced by £*. 

Remark 2.2.8. (i) Here in comparison with old estimates we replaced factor 
K2 by k which is an advantage; 

(ii) These estimates imply that S{ x -i(x)<\} l^l 2 dx < C/t 2 |log/?| while in 
fact it must not exceed Ck 2 . 

2.2.4 Estimates to a minimizer. Ill 

Consider now external zone y := {x : £(x) > 1} and assume that 

(2.2.32) CM < C£(x)~ v as £{x) > 1 
with v > 1. 

Then if also \dA(x)\ = 0{£{x)~^) as £{x) > 1 then the right hand 
expression of (2.2.29) does not exceed Ck{£^ v ^ + £^^- 2u ^ and therefore 
we almost upgrade estimate to |<9/4(x)| to 0{£~^ v + l^ u ^- 2v + 1 ^ an d repeating 
these arguments sufficiently many times to 0(£- 3u ). 
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However, there are obstacles to this conclusion: first, as v > 1 we 
conclude that 

Aj = J2*j,™\ x -y™r 1 + o(r 1 - s ) 

m 

with constant <Xj m ; however assumption V • A = implies <Xj m = and we 
pass this obstacle. The second obstacle 

A J = Yl Kjk.miXk ~ Yk,m)\>< ~ y m |~ 3 + 0(r 2 ) 
k,m 

with constant otjk, m we cannot pass as assumption V • A = implies only 
that modulo gradient A = J2 m |3 m x W" 1 with constant vectors |3 m and 
one cannot pass this obstacle. 

Therefore we upgrade (2.2.30)-(2.2.31) there: 

Proposition 2.2.9. In our framework assume additionally that (2.2.32) 
holds. Then as v > | 

(2.2.33) \A\ < CKt 2 , \dA\ < Ck£~ 3 
and 

(2.2.34) \dA{x) - 8A{y)\ < C e nr 3 - e \x - y\ e as \x - y\ < ^(x) 
as£{x) > 1 (for all 9 e {0,1)); 

Remark 2.2.10. (i) In application to the ground state energy we are inter- 
ested in v = 2; 

(ii) Observe that for a > 1 

(2.2.35) / \dA\ 2 dx= 0(/tV 3 ); 

(iii) We were not able to improve (2.2.33)-(2.2.35) no matter how fast ( 
decays. 
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2.3 Basic trace estimates 



Recall that the standard Tauberian theory results in the remainder estimate 
0(h~ 2 ). Really, as the effective magnetic field intensity is no more than 
Ck, contribution of B(x,£(x)) to the Tauberian error 13 ) does not exceed 
CC, 2 x hr 1 = C( 3 £h~ 1 which as ( x £^2 translates into C£~2h~ x and 
summation over {x : £(x) > £* = h 2 } results in Ch~ 2 . On the other hand, 
contribution of {x : £(x) < £* = h 2 } into asymptotics does not exceed 
Ch- 3 £~ l = Ch' 2 as h = 1. 

However now we can unleash arguments of V. Ivrii, V. and I. M. Sigal [IS]. 
Recall that we are looking at 

(2.3.1) Tr(V>H-^) = M<PiH AV <Pi) + Tr(0 2 H AV ,0 2 ) 

where 4> 2 = 4>\ + 02' su PP0i c {x, |x| < 2r}, supp0 2 C {x, r < \x\ < b} 
and we compare it with the same expression calculated for H A V o with 
V° = Z m |x| _1 . Here we assume that 

(2.3.2) a<l, zxl 
and 

(2.3.3) \D a (V - V°)\ < CqB-H-^ Va : \a\ < 3. 

The latter assumption is too restrictive and could be weaken. Then as 0(x) 
is an ^-admissible partition element 

(2.3.4) lr{Q{-H A V )<p 2 ) = J Weyl(x)0 2 (x) dx + 0(r/T 2 ) 

and 

(2.3.5) Jr(Hfy V (f) 2 ) = J \Ney\ 1 {x) ( f) 2 {x)dx+0{r-h- 1 ) 

where the error estimates are 0(h~ 2 ) and O^fir 1 ) respectively. One can 
justify transition from the Tauberian to Weyl errors by considering Tauberian 
expressions and considering Ha b ,v and Ha,v as unperturbed and perturbed 
operators respectively; their difference is 0(( 3 e 2 ) with e = h 1 ' 6 . 

Then the contribution 14 ) of time interval {t : t x T} to Tauberian 
expression for (2.3.4) of the first term in the approximation does not exceed 

13 ) Or Weyl error as we will explain transition from Tauberian to Weyl estimates below. 

14 ) After standard rescaling x i->- xl^ 1 , £ M> £C _1 , h ^ h, r i-> t(~ 2 , and t M> tC,(.~ x . 
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Ch~ A T x (hT^ 1 ) 5 , of the second term C7j~ 4 T x (hT~ 1 ) s Th~ 1 e 2 , and of the 
third term Ch~ A T x (hT~ 1 )T 2 h~ 2 e A . In the end the first tern gives us Weyl 
expression, the second term turns out to be 0, and the third term is less 
than the announced error. 

Similarly, the contribution 14 ) of time interval {t : t x T} to Tauberian 
expression for (2.3.5) of the first term in the approximation does not exceed 
C( 2 h~*T x {hT- 1 ) 5 , of the second term C( 2 h~ A T x (hT~ l ) 5 Th- l e 2 and of 
the third term Ch~ 4 T x (hT~ 1 ) 2 T 2 h~ 2 e 4 . Again, in the end the first tern 
gives us Weyl expression, the second term turns out to be 0, and the third 
term is less than the announced error. 

The same estimates also hold for operator Ha yo and then using £- 
admissible partition of unity we conclude that 

(2.3.6) Jr(MH AV -H Av0 ) ( j )2 ) = 

J (Weyl^x) - Weyl?(x)) <j>\(x) dx + O^h- 1 ) 

where Weyl? and Weyl are calculated for operator with potential V°. Indeed, 
we just proved this for each operator Ha,v and H A yo separately. 

On the other hand, considering V v = V°(l — rj) + Vrj = V° + Wrj 
and following V. Ivrii, V. and I. M. Sigal [IS] we can rewrite the similar 
expression albeit for 2 = 1 as 

(2.3.7) Tr / We{-H A ,vv)dr) 

Jo 

and applying the semiclassical approximation (under temporary assumption 
that W is supported in {x : |x| < 4r}) one can prove that for 0i = 1 

(2.3.8) Tr(0 1 (H- u -H- vO )0 1 ) = 

J (Weyl^x) - WeylS(x)) <j>\{x) dx + O^- 1 ^ 2 ). 

Really, due to (2.3.4) and (2.3.3) the contribution of ball B(x,£(x)) does not 
exceed Ca~ x ftr 2 = Ca~ 1 £(x)h~ 2 and summation with respect to partition as 
£(x) < 4r returns Ca~ 1 rh~ 2 ); meanwhile contribution of {x : £(x) < £*} does 
not exceed Ca~ l fT 2 = Ca^ 1 as there H=l. 

One can get easily rid off the temporary assumption and take 4>i sup- 
ported in {x : £{x) < 2r} instead. 

Therefore we arrive to 
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Proposition 2.3.1. Under assumption (2.3.3) 

(2.3.9) Tr(ij(H A V - H AV0 )^) = 

J (Weyl x (x) - Weyl?(x)) ^ 2 (x) dx + O^/r*) 

Really, a~3/j~5 is r^/j -1 + a~ x rbr 2 optimized by r x r, := (a/?)§; as 
h 2 < a we note that h 2 < r* < a. 

Corollary 2.3.2. (i) As M = 1 equality (2.3.9) remains valid with if) = 1 
and a = l. 

(mJ As M > 2 and a > h 2 equality (2.3.9) becomes 

(2.3.10) Tr(V(H- v - H^ty) = 

/ (Weyl^x) - Weyl?(x)) ^ 2 (x) dx + 0((a _ s + l)/rs) 

where we reset case a > 1 to a = 1. 

2.4 Improved trace estimates 
2.4.1 Improved Tauberian estimates 

Let us apply much more advanced arguments of Section 12.5 of [Ivrll]; 
recall that these arguments are using long term propagation of singularities. 
Unfortunately using these arguments we are not able to improve the above 
results unless k < 1. 

First, let us consider ip which is r-admissible partition element located 
in {x : £(x) x r} and we need to estimate an absolute value of 

(2.4.1) F t ^- lT xr(0Tr(e^ ltH ^) 

and to do it we need to estimate the same expression with Xr{t) replaced by 
XT'{t) with to < T' < T with t = <^C -1 = er5 - We can break ip = ip + + ip~ 
with ?/' ± = ^(x, hD) such that the trajectories in the positive (negative) 
time direction from support of its symbol i/j + (x, £) are going after time Qq 
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in the direction of increased £(x) and since we consider trace we need to 
consider only ip + and only x ^^([^ 1])- 

The trouble is that we have not rough but non-smooth magnetic field 15 ); 

3 

so let us consider to + ti + ... + t n x V where tj = er? , rj = c J r, j = 0, 1, ... , n 
and estimate an error appearing when we replace in (modified) (2.4.1) 

e ih- 1 tH q p+ by 

(2.4.2) e ^-tn)H^ e ih-H nH ^+ . . . e ih-h 1 H^+ e ih-%H^ + 

with tpf defined similarly and Hamiltonian flow from supp(^) for t = tj is 
inside of {(x, £) : ipj~ +1 (x, £) = 1}. Therefore we need to estimate an error 
when we insert 

According to our propagation results (namely, Proposition 1.2.2) after 

ipi , ... were inserted, insertion of ^ brings a relative error not ex- 

_ i 

ceeding C(^j , |||9/4||| 6) Y . + Hj + ) where hj = hr- 2 and Yj is an ery- vicinity of 
x-projection of supp(^); s is an arbitrarily large exponent. 

Recall that Y . < CkHj" 6 as 8 e (1,2); therefore this relative error 

does not exceed CHj(k + hj). So inserting all brings a relative error 
C^2j >0 hj(K + hj) x Ch(K + h s ) and since a priory expression (2.4.1) is 
bounded by Ch~ 3 T we conclude that 

(2.4.3) An absolute value of expression (2.4.1) with Xr{t) replaced by Xr>{t) 
with J"* x ri < V < T* 16 ) does not exceed OT 2 (k + h s )V '. 

Then an absolute value of expression (2.4.1) with Xr(t) replaced by Xt(0 — 
Xt, does not exceed Ch~ 2 (n + h s ) and since expression (2.4.1) with T = T* 
does not exceed Ch~ 2 t we conclude that 

(2.4.4) An absolute value of expression (2.4.1) with T* < T < T* does not 
exceed Ch~ 2 T* + CH~ 2 {k + h 5 )T . 

Then 

(2.4.5) An error when we replace Tr(0(— Ha ,v)if) by its Tauberian expression 
with "time" T does not exceed Ch~ 2 {T* T^ 1 + k + h s ) 

15) Q r roygh with the roughness parameter h which is a bit too short. 

16 ) We discuss the choice of T* later. 
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and 

(2.4.6) An error when we replace Tr( H^ v ifA by its Tauberian expression 
with "time" T does not exceed C/T 1 ( 7"* 7"" 1 + k + ft s ) 7"* T^ 2 . 

In the latter statement we need to remember how everything scales. 

Observe that presence of the magnetic field due to its estimates relatively 
perturbs dynamics by 0(k) and therefore as k is sufficiently small does not 
affect 7"*. Then assuming that 

(2.4.7) \V a {V- V°)\ < ea- l r~\ a \ Va : |a| < 1, V° = Zr 1 

with Z x 1, a > h 2 - s 

3 

we can take T* x a 2 and therefore we conclude that The Tauberian error 
in (2.4.6) does not exceed 

(2.4.8) C/j _1 a~5r(r5a _ i + k + h s r* s ) 

and we arrive to statement (i) below. 

Meanwhile The Tauberian error in Tr^(H^ v — H~ v0 )%Ij) does not exceed 

(2.4.9) Ca- l h- 2 r{Aa- Z 2 +K + h s r~^ s ) 
and we arrive to statement (ii) below: 

Proposition 2.4.1. Assume that (2.4.7) is fulfilled and let if) be I- admissible 
function supported in {x : £(x) x r} with h 2 < r < a. Let A satisfy 
minimizer estimate. Then 

(i) The Tauberian error with T = T* x at in Tr( H^ v if\ does not exceed 
(2.4.8); 

(ii) The Tauberian error with T = T* x a? in Tr((/-f^y — H^ v0 )ifj) does 
not exceed (2.4.9). 

Proof. An easy proof following arguments of Section 12.5 of [Ivrll] is left 
to the reader. □ 
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Note that summation in (i) with respect to r : b < r < a returns 
Ch~ 1 a~z and summation in (ii) with respect to r : h 2 < r < b returns 
Ch~ 2 (a~2 b? + Ka^b) + Ca~ x . Note also that Statement (ii) remains true 
for /--admissible function supported in {x : £(x) < r} with r x h 2 . Then we 
arrive to 

Corollary 2.4.2. Assume that (2.4.7) is fulfilled. Then 

(i) Let 02 be £- admissible function supported in {x : b < £(x) < a}. Then 
the Tauberian error with T = T* x as in Tr(H A v (f) 2 ) does not exceed 

Ch^a-h; 

(ii) Let 0! be £-admissible function supported in {x : £(x) < b}. Then the 
Tauberian error with T = T* x as in Tr((H AV — H AV o)(pi) does not exceed 
Ch- 2 {a--2bl + Ka^b) + Ca-\ 

Remark 2.4-3. Obviously we do not need any new assumptions on k to 
estimate the sum of expressions obtained in Statements (i) and (ii) of 
Corollary 2.4.2 by C/? _1 a _ 2 (as b < a^h) here but we need to move from 
Tauberian expression to Weyl' expression. 

2.4.2 Improved Weyl estimates 

Note that in virtue of (2.4.3) for element ip contribution of the time interval 
{t : \t\ x 7"'} to the Tauberian expression for Tr(^H A v tp} does not exceed 
Cfi _1 (/t + h s ) 7~* T'^ 1 ^ 2 and therefore replacing Xr(t) by Xr»(t) introduces 
an error not exceeding 

(2.4.10) CfT 1 (k + h s )( 2 x C/7-V-i (k + h 5 ^ 5 ) 
and summation with respect to r : b < r < a returns 

(2.4.11) Ch- l b- 1 *{K + h'b- 1 * 5 ). 

On the other hand, also in virtue of (2.4.3) for element ip contribution of 
the time interval {t : \t\ x 7"'} to the Tauberian expression for Tr((H AV — 
H AV0 )ip) does not exceed CtT 2 a~ l [k + h s ) and therefore replacing Xr(r) by 
XT,{t) introduces an error not exceeding Cbr 2 a~ x r [k + /rY~5 s )| log T/TJ. 
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Further, in virtue of (2.4.4) the Tauberian error does not exceed Ch 2 a 1 r[r^T 1 + 
k + /) s r _ 5 s ) and adding these two errors together optimizing the sum by 

3 5 1 i 

7 < 32 we get 7 x rz(/t + tfr'^y 1 and the sum 

(2.4.12) Ch- 2 a- 1 r{K + h s r- 1 2 5 ){\\og{K + h'r^ + l) + OrVM. 
Meanwhile repeating arguments of Subsection 2.3 one can see easily that 

(2.4.13) The difference between Tauberian expression with 7=7* and 
Weyl expression for Tr(H^ v tfj) does not exceed (2.4.10) with any s < 2 

and 

(2.4.14) The difference between Tauberian expression with 7=7* and 
Weyl expression for Tr((H^ u — H~ v0 )ip^J does not exceed 

(2.4.15) C/T 2 a- 1 r(fi; + /) s r~5 s ) 
with any s < 2 and thus does not exceed (2.4.12). 
Summation with respect to r : h 2 < r < b of (2.4.12) returns 

(2.4.16) Ch- 2 a~ x bK\ logre| + Ch^bh'^ + Ca -1 ; 

adding expression (2.4.11) and optimizing the sum by b : h 2 < b < a we get 

2 

b x (a/?| log/t|)3 and expression 

(2.4.17) C/7-t a -^| | g/t|3 + C/rV*. 
Thus we have proven 

Proposition 2.4.4. (%) In the framework of Proposition 2.3.1 

(2.4.18) Tr(^H^-H- u0 )V>) = 

y (Weyl^x) - Weyl°(x)) ^ 2 (x) cbc + 0(/|-3a-i«| log «|§ + r^-i). 
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(ii) In particular as 

(2.4.19) k < ca~h*\ loga/? 2 | 5 

the error in (2.4.18) does not exceed Ch~ 1 a~^ as in the case without magnetic 
field. 

Remark 2.4-5. (i) Obviously we could consider a = 1 and then just rescale 
x i — y xa -1 , r i — y ra, h i— >■ /?a~2; 

(ii) One may wonder if the same approach works for estimate of A. First 
of all, there is no improvement for estimate for \d 2 A\ as it follows from the 
estimate for \AA\ which is a pointwise estimate; 

(iii) However as dA and A are mollifications of AA one can improve estimates 
for them as k <C 1 and f <3; however there are no improvements as either 
k x 1 or t > a. Since these improvements do not lead to the improvements 
of our final results we do not pursue them. 



2.5 Single singularity 
2.5.1 Coulomb potential 

Consider now exactly Coulomb potential: V = Z|x| _1 . Let us establish the 
existence of the Scott correction: 

Proposition 2.5.1. Let V = Z|x| _1 ; h > 0, Z > and < k < k*. Then 

(i) The following limit exists 

(2.5.1) lim (inf(Tr(^H AV ^)" + ^ / \ dA \* dx ) 

-J Weyl^x^x) dx^j =: 2Z 2 /T 2 S(Z/t); 

(ii) And it coincides with 



-J \Ney\ 1 (H AV + r ] ,x)dx^ 
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(Hi) And with 

(2.5.3) inf^y (e 1 (H AV ;x,x,0)-\Ney\ 1 (H AV ,xf)dx+ 

P \dA\ 2 dx ) ; 



1 

^2 



(iv) We also can replace in (i) Tr(0 r /-//\_\/0 r ) by Tr(cf) r H A v (p r )- 
Here e ^ °°(e(0, 1)) ; = 1 in 6(0, |) ; r = 0(x/r). 

Proof. Due to scaling x h-> Zh~ 2 x, A h-> 2 r ~ 1 hA, dA h-> Z~ 2 h 3 dA one needs 
to consider only Z = /? = 1; all expressions on the left scale exactly as 
Z 2 /)- 2 S(Z/t). 

(i) Let us compare 

Q(r, re, /A) := Tr{<j> r H Ay <j> r )- - J \Ney\ 1 {H Ay , x)0 2 (x) + ^ |<9/l| 2 dx 
and Q(r', re, /4) with r > 1 and r' > 2r. Note that 

Q(r',K,A) > Q(r,(l + e)K,A)+ 

(M^2irHA,vM-- [ Weyl 1 (H AV ,x)^ v (x)dx+^- / |<94| 2 ^ 2 v (x) dx 
i<j<j "* K J 

where ip and ^ are smooth compactly supported functions, equal in 6(0, 2) 
and ip = 1 in the vicinity of supp(^), J = [\og 2 r'/VJ, e > is arbitrarily 
small. 

Therefore we can replace in the sum A by Aj and ^ 2 Jr by 1 in f \dA\ 2 ^> 2 jr (x) dx; 
but then in the virtue of Section 1 each term in the sum is bounded from 
below by -C'(e) / p 3 £~ 1 $ t dx = -C'(e)H with t = 2 j r. Then 

(2.5.4) Q(r', re, A) > Q(r, (1 + e)re, /A) - C'(e)H . 

We know that (?(r, (1 + e)re, /4) is bounded from below by — C(r) but now 
we conclude that this bound is uniform with respect to r. Then 2S(k) := 
lim inf r ^ +00 inf ^ Q(r, re, A) > —00. Further, (2.5.4) implies that 

2S(k) + C'(e)r-i > inf Q(r, (1 + e)re, 4) 
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and therefore 

(2.5.5) lim sup inf Q(r, k, A) < 2S((1 + e)" 1 ^). 

Furthermore, plugging A = we can see easily that Q(r, k, A) is uniformly 
bounded from above and therefore 2S(/t) < +00; also our arguments imply 
that J \dA\ 2 dx is uniformly bounded for near optimizers and therefore S(k) 
is continuous with respect to k < k*; combining with (2.5.5) we arrive to 
Statement (i). 

(ii) Similarly, (2.5.4) holds with H AV replaced by H AV + r\ and then we can 
take r' = 00 and apply inf ^ to both sides arriving to 

\nj(Tr(H A y + riy - J Weyl^,^, x) dx + ^ J \dA\ 2 dx^j > 

inf(Vr(0 r (H AU + T7)0 r ;r -J \Ne ) /\ 1 {H A y +r) , x)^(x) dx+ 

jz 1 , / \dAAdx-C{e)r*. 
(1 + e)Kj J W 

After this as rj — > +0 the right-hand expression tends to itself with 77 = 0; 
tending r — > +00 we get there 2S((1 + e)/t) in virtue of Statement (i) and 
tending e — > +0 we arrive to 

(2.5.6) lim inf inf (jr(H AV + ri)~ - J Weyl^H*, V+r/, x) dx+ 

^ j \dA\ 2 dx^j > 2S(k). 

On the other hand, consider 

Tr(H AV + v y - J \Ney\ 1 (H AV+v ,x)dx+^ J \dA\ 2 dx 

and replace Tr(H A y + rj) by 

(2.5.7) Tr((H Ay + r,y<t> 2 ) + -^ J \dA\ 2 dx + Tr[[H A y + ~ 0r)) ■ 
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Let A be a minimizer of the first expression; then in virtue of Propo- 
sitions 1.3.4-L3.6 this minimizer is sufficiently "good" on er-vicinity of 
supp(l — <p 2 r ) that the the difference between the second term and its Weyl 
expression does not exceed Cr~2; one can prove it easily by £(x)-admissible 
partition of unity as in part (i) of the proof and we leave details to the 
reader. 

Observe that the first term in (2.5.7) is inf^Tr^^H^y + 77) r )- In 
this expression we can take limit as rj — > +0 just setting 77 = and therefore 
the left-hand expression in (2.5.6) with lim inf replaced by limsup does not 
exceed 

inf(Tr(0 r H^ r ) - y Weyl^H^, x)<f> 2 (x) dx + j \dA\ 2 dx) + Cr~h 

taking limit as r -> +00 we conclude that the left-hand expression in (2.5.6) 
with lim inf replaced by limsup does not exceed 

liminfinffTr(0 f H7w0 r ) - / Weyl^H* v, x)(f) 2 {x) dx + - [ \dA\ 2 dx) . 

A \ ' J K J J 

This expression does not exceed (2.5.1) and therefore combining with (2.5.6) 
we prove Statements (ii) and (iv). 

(iii) Similarly 

J (e 1 (H A u;x,x,0)-Weyl 1 (H A v,x))^(x)dx+i J \dA\ 2 dx > 

mf(Tr(<f> r H AiV (l> r )~ - J Weyl 1 (x)^(x) dx + J \dA\ 2 ^ - Cr* 
and therefore 

inf lim inf / (e x {H A v ; x, x, 0) - Weyl^H^ v , x))<f> 2 r {x) dx+ 

A r^oo J \ ' J 

^ J \8A\ 2 dx > 2S(k). 

On the other hand, as in (ii) taking A to be a minimizer of the first expression 
in (2.5.7) we see that 

inf limsup / (e 1 (H AV ;x,x,0)-\Ney\ 1 (H A y,x))(p 2 (x)dx+ 

^ y \dA\ 2 dx < 2S{k) 
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and Statement (iii) has been proven. □ 

Remark 2.5.2. (i) Statements similar to(i), (ii) were proven in L. Erdos, 
S. Fournais, and J. P. Solovej [EFS3] (see Theorem 2.4 and Lemma 2.5 
respectively). 

(ii) Again as observed in in L. Erdos, S. Fournais, and J. P. Solovej [EFS3] 
we do not know if (a) S{k) < S(0) for k > or just (b) S(k) = S(0) as 
k < k* and 5(/t) = — oo as k > k*. If we knew that the optimizer is unique 
then obviously A = and it would be relatively easy but rather unexciting 
the latter case. 

(iii) While we assumed that k < k* with k* > and it is possible that 
S(k) = — oo as k > k* with some k* < oo we are not aware about any proof 
of this, so in fact it could happen that k* = +oo and then condition k < k* 
is superficial and one needs to study asymptotics of S(k) as k — )■ +oo. 

Proposition 2.5.3. As < k < k' 

(2.5.8) S{k') < S{k) < S{k') + Ck\k- 1 - k'- 1 ). 

Proof. Monotonicity of S(k) is obvious. 

Let < k < k,' < k" < k* . Then for any e > if r = r £ is large enough 
then the left-hand expression in (2.5.1) for k' (without inf and lim) is greater 
than S(k") — e + (k'^ 1 — k //_1 )||c}/4|| 2 ; also, if A is an almost minimizer there, 
it is less than S(k') + e. 

Therefore {k'^ 1 - K"- 1 )\\dA\\ 2 < \S{k") - S{k')\ + 2e. But then 

5{k) -e< 5{k') + e + (/sT 1 - « /_1 )||^|| 2 < 

5{k') + e+ C(/€~ 1 - k'- 1 )^'- 1 - k"- 1 )- 1 ^^") - S(k')| + 2e) 
and therefore 

(2.5.9) {k- 1 - k'- 1 )- 1 ^) - S«)| < {k'- 1 - k"~ 1 Y 1 \S{k') - S(k")| 
which for k" = k* implies (2.5.8). □ 
Remark 2.5.4- Using global equation (1.1.14) we conclude that for Z = h = 1 

(2.5.10) \d a A\ < CKi' 1 -^ 1 as t > 1, \a\ < 1, 

(2.5.11) \d a A\ <Ck as t < 1, |a| < 1, 
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(2.5.12) \\dA\\ 2 < Ck 2 . 
Then 

(2.5.13) S'(k)<C, |S(k(1 + 77)) - 5{r])\ < Ckt). 

2.5.2 Main theorem 

In the "atomic" case M = 1 we arrive instantly to 
Theorem 2.5.5. Let M = 1 and k < k* . Then 
(i) Asymptotics holds 

(2.5.14) E* = J\Ney\ 1 {x)dx + 2z 2 S{zK)h' 2 + 0(Ai-i«j log/c|i + /T 1 ); 
(m,) = o(/?3| log /?| 3) £/ien 

(2.5.15) E* = J Weyl*(x) cbc + 2z 2 S(zK)h~ 2 + o^ 1 ). 

Proof. If /4 satisfies minimizer properties then in virtue of Proposition 2.4.4 

(2.5.16) Tr~ H A „ - y Weyl^x) dx = Tr~ H AV o - ^ Weyl°(x) cfx 

+ 0(/i~3k| log k| 3 + /T 1 ) 

and adding magnetic energy and plugging either minimizer for V or for V° 

we get 



(2.5.17) inf (Tr - H Ay - J Weyl^x) dx + J \dA\ 2 dx^j 
inf (Tr _ H AV o - ^ Weyl?(x) dx + J \dA\ 2 dx 



< 
> 



± C(/T3/c| log /t|3 + /T 1 ). 

Sure as V (and surely V°) are not sufficiently fast decaying at infinity the 
left (and for sure the right hand) expression in (2.5.16) should be regularized 
as in Subsection 2.5.1. However for potential decaying fast enough (faster 
than |x|~ 2 ~ 5 ) regularization is not needed. 
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For V° we have an exact expression which concludes the proof of State- 
ment (i). 

The proof of (ii) is similar albeit with the small improvement based on 
the behavior of the classical dynamics (without magnetic field) exactly as in 
Chapter 24 of [Ivr 11]. □ 

2.6 Several singularities 

Consider now "molecular" case M > 1. Then we need more delicate 
arguments. 

2.6.1 Decoupling of singularities 

Consider partition of unity 1 = X]o<m</w^m where ip m is supported in la- 
vicinity of y m as m = 1 , . . . , M and ipo = in t a- vicinities of y m ( "near- nuclei" 
and "between- nuclei" partition elements). 

Estimate from above. 

Then 

(2.6.1) Tr H^ v = M^H- v^) 

0<m<M 

and to estimate E* from the above we impose an extra condition to A: 

(2.6.2) A = as £(x) > -a. 

5 

Then in this framework we estimate 

(2.6.3) Tr(^oW^o) - J WeyliWoM dx < Ch^a^. 

Proof of this inequality is trivial by using ^-admissible partition and applying 
results of the theory without any magnetic field. 

So, to estimate E* from above 17 ) we just need to estimate from above 
the minimum with respect to A satisfying (2.6.2) of the expression 

(2.6.4) Jr{i, m H AV i, m )- /weyl 1 (x)^(x)c/x+4^ / \dA\ 2 dx. 
17) Modulo error in (2.3.10). 
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Estimate from below. 

In this case we use the same partition of unity {^mli^o.i m an d estimate 

(2.6.5) TV H A v > ^ Jr{^ m H AV ^ m ) 

0<m<M 

with 

(2.6.6) V = V + 2h 2 (^) 2 

0<m<M 

and we also use decomposition 

(2.6.7) j\dA\ 2 dx = [ulldAfdx 
with 

(2.6.8) u m {x) = 1 as £ m {x) < -^a, u m {x) >1-Cq as £ m {x) < l -a 

m = 1 M, 

(2.6.9) w > e ^ as £(x) > ^a. 

5 

So far <r > is a constant but later it will be a small parameter. Then since 

(2.6.10) Tr-(^oH A v"^o) " J Weyl 1 (x)^(x) dx + ^ J " 2 \dA\ 2 dx > 

(again proven by partition) in virtue of the previous Section 1 we are left 
with the estimates from below for 

(2.6.11) Jr(^ m H A , v ^ m ) - J Weyl 1 (x)^(x) dx + J u 2 m \dA\ 2 dx. 

Remark 2.6.1. (i) Note that the error in Wey^ when we replace V there 
by V does not exceed C/j _1 (1 + a^s) which is less than the error in (2.3.10). 
Here we can also assume that A satisfies (2.6.2); we need just to replace q 
by e <? in (2.6.8)-(2.6.9). 

(ii) We can further go down by replacing Tr~ (ip m HA,v>ip m ) by ^f{ipmH A v ,ip m )- 
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(iii) Therefore we basically have the same object for both estimates albeit 
with marginally different potentials ( V in the estimate from above and V in 
the estimate from below) and with a weight u 2 m satisfying (2.6.8)-(2.6.9); in 
both cases u = 1 as £(x) < ^a but in the estimate from above cu(x) grows 
to C and in the estimate from below u(x) decays to q as i(x) > |a and in 
both cases condition (2.6.2) could be either imposed or skipped. 

(iv) From now on we consider a single singularity at and we skip subscript 
m. However if there was a single singularity from the beginning, all arguments 
of this and forthcoming paragraphs would be unnecessary. 

Scaling. 

(i) We are done as Z x 1 but as Z < 1 18 ) we need a bit more fixing. 
The problem is that V x Zi^ 1 only as |x| < aZ; otherwise V < a -1 
(where we still assume that a < 1). To deal with this we apply in the zone 
{x : aZ < |x| < a} the same procedure as before and its contribution to 
the error will be Ch~ l a~z as p = a~2 here. Actually we also need to keep 
|x| > Z -1 /) 2 ; so we assume that Z -1 /? 2 < Za i.e. Z > a - 2/7. 

Now let us scale x h-» x' = x(aZ) -1 , and multiply H a V by a and therefore 
also multiply A by aa, so A i-> A' = aM, h h-> h' = ha^iZ^ 1 ; then 
the magnetic energy becomes /t~ 1 /?~ 2 Z J u 2 (x)\d' A'\ 2 dx' where factors a -1 
and aZ come from substitution A = a~iA' and scaling respectively. We 
need to multiply it by a (as we multiplied an operator); then plugging 
h~ 2 = h'~ 2 a~ 1 Z~ 2 we get the same expression as before but with Z' = 1, 
a' = 1 and h' = ha~zZ~ l < 1 and k' = kZ instead of h and k. 

If we establish here an error 0(/7 / - 1 + /t / |log/t / |3/ 7 / -3) the final error will 

be this expression multiplied by a" 1 i.e. 0{aT^Zbr x + na~^Z^ | log/tZ|3 h~i) 
which is less than the same expression with Z = 1. 

(ii) On the other hand, let Z < a~^h. Recall, we assume that a > Coh 2 . 
Then we can apply the same arguments as before but with Z = a~^h and we 
arrive to the same situation as before albeit with h' = 1, a' = 1, k! = na~ih 
and with Z' = ZZ l . Then we have the trivial error estimate 0(a _1 ) which 
is less than 0(a^2h^ 1 ). 

18 ) As Z denotes Z m now we assume only that Z\ + ... + Zm x 1. 
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2.6.2 Main results 

Combining results of the previous Subsubsections and Paragraphs with 
proposition 2.2.9 we arrive to 

Theorem 2.6.2. Let M > 2, k < k* and (2.2.32) hold, with v > |. Then 
(i) Asymptotics holds 

(2.6.12) E* = f Weyl^x) dx + 2 ^ Z 2 m S{Z m K)h- 2 + 0(/?! + R 2 ) 

J l<m<M 

with 



(2.6.13) /?i = 
and 

(2.6.14) /? 2 = nh~ 2 



_1 1 4 

/? + /c| log k| 3/7 3 as a > 1 

— - — 1 ' 1 - — - — - 2 
32/7 + K | log K\ia 3/) 3 as /? < a < 1 

1 

a as a > \ log h\ 3, 

I log/iV 1 ! -1 as /? 2 < a < I log/»|3; 



fnj If k, = o(/?3 1 log /j| 3) ; Ka 3 = o(/?) and a 1 = o(l) inen 

(2.6.15) E* = / Weyi;(x) dx + 2 Z 2 m S{Z m K)h' 2 + o(/ ? - 1 ). 

l<m<M 

Proof. To prove theorem we need to prove the following estimate 

(2.6.16) hdA\\ 2 {b<f{x)<2b} < Cb- 3 

where r* < b < a is a "cut-off" . On the other hand we know that 

1 (99 

(2.6.17) -||«9/l|| 2 = -— = 0(1) 

and we need to recover the last factor in the definition of R 2 . 

As a > 1 we can have na~ z because in virtue of (2.2.35) the square of 
the partial norm in (2.6.17) does not exceed Ca~ 3 n 2 . 

On the other hand, as h 2 < < a we can select b : r* < b < a such that 
he partial norm in (2.6.17) does not exceed C| log(a/ /7 2 )| 1 • ||<9/4|| 2 . □ 
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Remark 2.6.3. (i) As a < | \ogh\ we do not need assumption (2.2.32); 

(ii) Following arguments of Section 12.5 of [Ivrll] estimates (2.5.15) and 
(2.6.17) could be improved to 0(h~ 1+s ) provided a > h~ Sl , k < hi +8l \ log /7| 5 
and k < a 3 h 1+Sl . 

2.6.3 Problems and remarks 

Problem 2.6.4. (i) As k G [0, k*] with small enough k* > does S(k) 
really depend on k or S(k) = S(0) (see Remark 2.5.2)? 

(ii) If S(k) really depends on k, what is asymptotic behavior of S(/t) — S(0) 
as K — y +0: can one improve S(k) — S(0) = 0(k) ? 

(Hi) Do we really need an assumption k G [0, k*] (again see Remark 2.5.2)? 
(iv) Can one improve estimates to minimizer? 

3 Asymptotics of the ground state energy 
3.1 Problem 

Now we are ready to tackle our original object (0.0.2)-(0.0.3). So, let us 
consider our usual quantum Hamiltonian 



(3.1.1) 


H= E "? + 


Yl \ X J- X k\ 1 




l<j<N 


l<j<k<N 


in 






(3.1.2) 


*= A H ' 

l<n<W 


H = ^ 2 (M 3 ,C 2 ) 


with 






(3.1.3) 


H° = ((/V - 


A) -a) 2 - V(x) 



We are interested in the ground state energy E^(/4) of our system i.e. in the 
lowest eigenvalue of the operator H on S): 

(3.1.4) E* w (0) = inf Spec(H) on Sj 



59 



as A = and more generally in 

(3.1.5) = inf (mf Spec fl (H) + i J |V x 4| 2 cbc) 
where 

(3.1.6) V(x) Yl Zm 



\x-y 



l<m<M 

(3.1.7) WxZ»l, Z := Zi + ... + Z M , Z x > Z M > 

M is fixed, under assumption 

(3.1.8) < a < ^Z" 1 

with sufficiently small constant > 0. 
Our purpose is to derive an asymptotics 

(3.1.9) e;^4 f + 2Z l s ^) 

l<m<M 

and estimate an error (usually) provided 

(3.1.10) b := min |y m - y m ,| > Z~l 

l<m<m'<M 

Recall that Thomas-Fermi potential 1/1/ TF and Thomas-Fermi density 
p TF satisfy equations 

(3.1.11) p TF = ^(l/l/ TF + v)\ 

and 

(3.1.12) l/1/ TF = V°+ |x|- 1 *p TF 
where v is a chemical potential. 
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3.2 Lower estimate 

Consider corresponding to H quadratic form exactly as in Sections 24.2 
and 25.6 of [Ivrll] 

(3.2.1) <Hv|/,v|/)=^(tf°v|/,v|/) + ( \xj-x k \- 1 V,V) = 

j l<j<k<N 
j l<j<k<N 

with 

(3.2.2) H = ((/'V - A) • tr) 2 - l/l/(x) 

where we select 1/1/ later. By Lieb-Oxford inequality the last term is estimated 
from below: 



(3.2.3) ( k -Xfcl^vi/.M/) > iD(p WlPv )- C /"p|c/x 

l<j<k<N J 

where 

(3.2.4) py{x) = N J \^{x,x 2 ,...,x N )\ 2 dx 2 ---dx N 
is a spatial density associated with v|/ and 

(3.2.5) D(p, p') ■= Jj\x- Y l~ V(x)p'(y) c/xc/y. 
Therefore again repeating arguments of Section 24.2 of [Ivrll] 

(3.2.6) <Hvj/,vj/)> 

^(H x .v|/,v|/)-2((V- W/)v|/,v|/) + iD( Py „Pvu)-C j p\dx = 
j 

^(H X; M/, M/) - D(p, p„) + ^D(pv, p*)-C f pi dx = 
j J 

^(H Xj M/, M/) - ^D(p, p) + l -D(p - pv, p-pv)-C j p\dx 
j J 
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as 

(3.2.7) W-V = \x\~ 1 *p. 
Note that due to antisymmetricity of U' 

(3.2.8) E( H ^- > E \->Tr"(H) 

J l<j<N:Xj<0 

where Ay are eigenvalues of /-/. 

To estimate the last term in (3.2.6) we reproduce the proof of Lemma 
4.3 from L. Erd6s, S. Fournais and J. P. Solovej [EFS3]: 

According to magnetic Lieb-Thirring inequality for U > 0: 

(3.2.9) E - U)^. V)>-C [ U 5/2 dx - C-f- 3 U A dx - 7 I B 2 dx 

j<N J J 

where B = V x 4, 7>0is arbitrary. Selecting U = f3 min(p^ / //3 , 7P^/ 3 ) with 
f3 > small but independent from 7 we ensure \Up^ > CU 5 ^ 2 + C"f~ 3 U 4 
and then 

(3.2.10) E <( H ?-)V, ^)> e [ min (pS/ 3 . lp A/3 )dx - 7 j B 2 dx 

j<N J J 

which implies 

(3.2.11) J p^ 3 dx < 7- 1 J m\u{p 5 J 3 , 1P ^)dx + 7 J p^dx < 

C7" 1 E (( H xj)V, V)+c f B 2 dx + c 7 /V 

j:Xj<0 J 

where we use J pydx = N. 

Remark 3.2.1. As one can prove easily (see also L. Erdos, S. Fournais and 
J. P. Solovej [EFS3]) that 

(3.2.12) E (i H i)^' V) < CZ ^ N 

j<N 
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even if N ^ Z; then we conclude that 

(3.2.13) J p 4 J 3 dx < CZi/V + Q / B 2 dx. 

It is sufficient unless we want to recover Dirac-Schwinger terms which 
unfortunately is possible only as a C l°g Z\~i. To recover remainder 
estimate o(Zs) (or marginally better) we just apply Theorem 25. A. 2 of 
[Ivrll]. We will do it later (see Theorem 3.1.5). 

Therefore skipping the non-negative third term in the right-hand expres- 
sion of (3.2.6) we conclude that 

(3.2.14) (HV, V) + - / |V x A\ 2 dx > 

a J 

Tr{H) + (- - d) / |V x A\ 2 dx- lD(p,p) - czl 
a J 2 

Applying Theorem 2.6.2 we conclude that 

(3.2.15) The sum of the first and the second terms in the right-hand expres- 
sion of (3.2.14) is greater than 

£ TF + ^2Z2S(aZ m ) - CZt(/? x + R 2 ) 

m 

with /?i and /?2 defined by (2.6.13) and (2.6.14) respectively with k = aZ, 
h = Z~3 and 



(3.2.16) a := Zs min |y m — y 

l<m<m'<M 



To prove this claim one needs just to rescale 

(3.2.17) ^5 x^xZJ, a^aZi W^Z-^W, 

A A, Vx/lhf ZsV x A 

and introduce 

(3.2.18) h = Z~K K = aZ. 
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4 

Observe that due to (3.2.17) 3 we need to multiply our estimate by Zs. 

Here one definitely needs the regularity properties like in Section 2 but 
we have them as p — p TF , 1/1/ = 1/1/ TF . Also one can see easily that "— Q" 
brings correction not exceeding C2«Z 2 as aZ < 1. 

Meanwhile for p = p TF , 1/1/ = 1/1/ TF 

(3.2.19) -2-Jw'*dx- ^D(p, p) = S TF . 

Lower estimate of Theorem 3.4.1 below has been proven. 

Remark 3.2.2. p = p TF , W = W TF delivers the maximum of the right-hand 
expression of (3.2.19) among p, W satisfying (3.2.7). 

3.3 Upper Estimate 

Upper estimate is easy. Plugging as in Section 24.2 of [Ivrll] the Slater 
determinant (24.2.16) of ip\, ... ,ipN where ... , ■?/>« are eigenf unctions of 
Ha,w we get 

(3.3.1) (HV|/, v|/) = Tr {H AW - X N ) + X N N+ 

J (1/1/ - V)(x)pv(x) dx + -D(pm,, pv)- 

^N{N-1) J J |xi -x 2 |- 1 |M/(x 1 ,x 2 ,x 3 x w | 2 dx x ■ ■ ■ dx N 

where we do not care about last term as we drop it (again as long as we 
cannot get sharp enough estimate) and the first term in the second line is 
in fact — D(p,pv) provided (3.2.7) holds. Thus we get 

(3.3.2) Tr(H AW - A w ) + X N N - ^D(p, p) + ^D(p„ - p, p v - p)+ 



1 

a 



\dA\ 2 dx 



where we added magnetic energy. Definitely we have several problems here: 
A/v depends on A and there may be less than N negative eigenvalues. 

However in the latter case we can obviously replace N by the lesser 
number N' := max(n < N, X n < 0) as is decreasing function of N. In this 
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case the first term in (3.3.2) would be just Tr~ (Ha ,w) and the second would 
be 0. Then we apply theory of the previous Section 2 immediately without 
extra complications. 

Consider A a minimizer (or its mollification) for operator Ha t w ~ P with 
potential 1/1/ = 1/1/ TF and p < 0. Then 

(3.3.3) N(/x) := #{A, < p) = J (W + p)\ dx + 0(Z§). 

One can prove (3.3.3) easily using the regularity properties of A established 
in the previous Section 2 and the same rescaling (3.2.17)-(3.2.18) as before. 
We leave this easy proof to the reader. 

Then repeating arguments of Subsubsection 24.4.2.1 "Estimating \\n— v|" 
of [Ivrll] we conclude that either N > Z — C Zi and then \v\ < CiZs and 
we can take p = and |A/v| < CiZs or N < Z — CqZ? and then we take 

i l 2 

H = v, | A a/ | x |z/| x (Z — A/)| and | A/v — ^| < Ci\u\*Z3. 

Then following again to arguments of Subsection 24.4.2 of [Ivrll] we 
conclude that 

(3.3.4) Expression (3.3.2) without term |D(py/ — p, — p) does not exceed 

£ TF + 2Z ^(«Z m ) + CZl(R, + /? 2 ) 

/T) 

with /?i and /?2 defined by (2.6.13) and (2.6.14) respectively with k = aZ, 
h = Z~5 and a defined by (3.2.16). 

Now we need to estimate properly D(py — p, py — p) which as in i.e. 
Subsubsection 24.4.2.2 "Estimating D-term" of [Ivrll] does not exceed the 
sum of 

(3.3.5) D(e(x, x, p) - Weyl(x, p), e(x, x, p) - Weyl(x, p)), 

(3.3.6) D(e(x, x, X N ) - Weyl(x, X N ), e(x, x, X N ) - Weyl(x, X N )), 

and 

(3.3.7) D(Weyl(x, p) - Weyl(x, X N ), Weyl(x, p) - Weyl(x, X N )). 

Following arguments of Subsubsection 24.4.2.2 "Estimating D-term" of 
[Ivrll] one can prove easily that due to regularity properties of A both two 

5 

semiclassical terms do not exceed CZs and due to estimates to | A/v — p\ the 

5 

last term does not exceed CZs either. 

This concludes the proof of the upper estimate in Theorem 3.4.1 below. 
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3.4 Main theorems 

Theorem 3.4.1. (i) Under assumptions (3.1.7) and (3.1.8) 

(3.4.1) E; = 4 F + 2Z 2 m S(aZ m ) + 0(zl(R 1 + R 2 )) 

l<m<M 

with /?i and R 2 defined by (2.6.13) and (2.6.14) respectively with k = aZ, 
h = Z~3 and a defined by (3.2.16), a = oo as M = 1; 

(ii) In particular under assumption (3.1.10) the following estimate holds 

(3.4.2) E* W = 4 F + Yl 2Z 2 m S(aZ m )+ 

l<m<M 

0(a\ log(aZ)|izf + Z§ + ab^Z 2 ) 

Recall that £jf is a Thomas-Fermi energy and S(Z m )Z 2 are magnetic 
Scott correction terms. 

Theorem 3.4.2. (i) Let assumptions (3.1.7) and (3.1.8) be fulfilled and let 
— v|/ A be a ground state for a near optimizer A of the original multiparticle 
problem. Then 

(3.4.3) D(p„-p TF , p„-p TF )<CZl; 

(ii) Furthermore, as b > Z~? 

(3.4.4) D(p„-p TF , Pw -p TF )< CZl(Z- s + (bZ 1 s)- 5 + (aZ) 5 ). 

Proof, (i) Note that all the terms in estimates from below and from above 
are OyZs) except the common term 

(3.4.5) Tr(H A w + aO + - [ |V x A\ 2 dx 

a J 

where A is a minimizer for this term and therefore estimate (3.4.3) has been 
proven because estimate from below also contains D(p^ — p TF , p^ — p TF ). 

(ii) To prove Statement (ii) one needs 
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(a) To improve estimate (3.3.3) to 

(3.4.6) N(/i) = J(W + fi)l dx+ 0{zl[Z- 5 + {bZ 1 i)- 5 + {aZf]), 

(b) To estimate terms (3.3.5)-(3.3.7) by the right-hand expression of (24.4.44) 
of [Ivrll], and 

(c) To accommodate Dirac term in both upper and lower estimates. 

Tasks (a), (b) are easy and we leave it to the reader (cf. arguments of 
Subsection 24.4.3) of [Ivrll]; we use that after rescaling effective magnetic 
field intensity becomes 0(aZ) in zone {x : £(x) x Z~5} due to already 
established estimates to A 

To fulfill (c) note that in the upper estimate we already have term 

(3.4.7) --tr J J \x-y\- 1 e N {x,y)e ] N {x,y)dxdy. 

On the other hand, in virtue of Theorem 25. A. 2 of [Ivrll] we replace in the 

4 5 

lower estimate term — Cf p^,(x) dx by (3.1.6) with 0(Z3~ 5 ) error (again we 
leave easy details to the reader). 

One can prove by the same arguments as as in the non-magnetic case 
that for aZ < k* it is Dirac + O(ZH). □ 

Finally, combining arguments sketched in the proof of Theorem 3.4.2 
with the improved estimate of (3.4.5) (see Theorem 2.6.2(h)) we arrive to 

Theorem 3.4.3. Let assumptions (3.1.8) and (3.1.10) be fulfilled, and let 
a < Z~t | log Z|~3 . Then 

(3.4.8) E* W = 4 F + J2 2Z^S{aZ m ) + Dirac + Schwinger+ 

l<m<M 

0(a\ log(aZ)|3Zf + zM + a b~ 3 Z 2 ) 

where Dirac and Schwinger are Dirac and Schwinger correction terms de- 
fined exactly as in non-magnetic case by (24.1.29) and (24.1.30) of [Ivrll] 
respectively. 
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3.5 Free nuclei model 

Consider now free nuclei model (see Subsubsection 24.4.4.2) of [Ivrll]. 
Theorem 3.5.1. Let us consider y m = y* m minimizing the full energy 

(3.5.1) E* N :=E* N + ^ Z m Z m ,\y m - y m '\~ 1 . 

l<m<m'<M 

Assume that 

(3.5.2) Z m x/V Vm=l M. 

Then 

(3.5.3) b > m\n(Z-^ +s , Z-^(aZ)- 6 , a^Z^) 

and in the remainder estimates in (3.4.2) and (3.4.8) one can skip b- 
connected terms; so we arrive to 

(3.5.4) E* W = 4 F + 2Z 2 m S(aZ m ) + 0(a\\og(aZ)\ 1 iZ 2 i + Z§) 

l<m<M 

and 

(3.5.5) E* N = S^ F + 2Z m 5 ( aZ ™) + Dirac + Schwinger+ 

l<m<M 

0(a|log(aZ)|3ZT +zH) 

respectively and also the same asymptotics with E* N and £jf instead of E* N 
and £jf . 

Proof. Optimization with respect to ... ,y M implies 

(3.5.6) E* + J2 UT^V zZ E ™ 

l<m<m'<M '" m " m '' l<m<M 

where E* = E*(y 1 , ... ,y M ; Z 1: ... , Z m , N) and E^ = E*(y m , Z m ) are calculated 
for separate atoms. In virtue of theorem 3.4.1 

(3.5.7) £ TF + £ r^^- £ ^ 

l<m<m'<M ' m m '' l<m<M 

Ca| log(aZ)|szf + Z§ + Cab^Z 2 ; 
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however due to strong non-binding theorem in Thomas-Fermi theory the 
left-hand expression is > b~ 7 as b > Z 3 and therefore (3.5.7) implies 

b > m\r\(Z~^, a -?| \og(aZ)\~^-Z~^ , aT*Z~%) 

where the third expression is larger than the second one for sure. Unfortu- 
nately for a > Z~~9~ s it is not as good as we claimed in (3.5.3). Still this 
estimate implies both (3.5.4) and (3.5.5). 

To prove (3.5.3) we record that b 3> Z ? 19 ^ we employ arguments used 
in the proof of Proposition 25.8.12 of [Ivrll] and prove that 

\Tr(H AW + 11) - ^ Tr ~i H A,w m 

l<m<M 

[ (Wey\{H A , w +[i; x)- ^ \Ney\{H A , Wm +n; x)) dx\ < CZ§ (Z- 5 +(aZ) 5 ) 

J l<m<M 

where A be a minimizer for "molecular" expression (3.4.7) and W m are 
atomic potentials. The same estimate holds if we replace Tr~ (Ha ,w m + by 
Tr~{H Am ,w m + H) with A m = A(j>(b~ x \x - y m |) with <\> e %°°{B{0^), equal 1 
in 6(0, |). We leave an easy proof to the reader. 

Then using lower estimate for inf Spec f:i (H) and upper estimates for 
inf SpeCj5 m (H m ) through Tr~ (Ha,w + /-O and Tr~(HA m ,w m + A*') respectively 
(where H m are associated with HA m ,w m ) we arrive to 

infSpe Cj5 (H)> infSpec^(H m ) + ^ TF - £ S T J 

l<m<M l<m<M 

- CZl (Z- & + (qZ) 5 ) 

and therefore 

(3.5.8) E A > E mAn + £ TF - 

l<m<M l<m<M 

CZl (Z- s + (aZ) s ) - Cab- 3 Z 2 

where the last term is due to replacement of - J |V x A\ 2 dx by "atomized" 
expressions Xli< m <M a 

J\V x A m \ 2 dx. 

19 ) There is no binding with b < Z~i because remainder estimate is (better than) CZ 2 
and binding energy excess is x Zs. 
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The last inequality (3.5.8) then obviously holds with A m replaced by 
optimizers for "atomic" expressions (3.4.7) and now strong non-binding 
theorem implies that 

b~ 7 < CZl(Z- 5 + (aZ) 5 ) + Cab- 3 Z 2 

which implies (3.5.3) where we change 5 > as needed. □ 



4 Miscellaneous problems 

In our analysis in Sections 24.5 and 24.6 of [Ivrll] the crucial role was played 
by an estimate of D(p^ — p TF , p-y — p TF ) and since what we have now (see 
Theorem 3.4.2) is (almost) as good as we had then, all arguments of these 
Sections still work with the minimal modifications. We leave most of the 
easy details to the reader but we need to deal with different magnetic fields 
for different N. 

4.1 Excessive negative charge 

Theorem 4.1.1. 20 ) Let condition (3.5.2) be fulfilled. 

(i) In the framework of the fixed nuclei model let us assume that 
'w : = E *a/-i " E* w > 0. Then 

5 f 1 as a < Z 

4.1.1 (N-Z)+<CZU . ~ . 

1 \ Z- s + (aZi)- s + (aZ) 5 as a > Z-l, 

(ii) In particular for a single atom and for molecule with a > Z~3+ s 
(4.1.2) {N-Z) + <Z 5 i(Z- s + {aZ) s ); 

(Hi) In the framework of the free nuclei model let us assume that \* N := 
E* N _! — Eli > 0. Then estimate (4.1.2) holds. 

2 °) cf. Theorem 24.5.2 of [Ivrll]. 
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Proof. The proof of follows the proof of Theorem 24.5.2 of [Ivrll] and since 

3 

it is not non-magnetic field specific we find that (N — Z)+ < where Q is 
an estimate for D(pvj/ — p TF , pm — p TF ) which we established already; recall 
that p TF x £~ 6 as £ > Z~3 also plays important role. 

Here we pick up A = A^ (exactly as in the analysis of free nuclei model 
we pick up y for A/-electrons) and conclude that 

\ N (A) := E N ^(A) - E N (A) > E* N _, - E* N > 

and then repeat arguments of the proof of Theorem 24.5.2 of [Ivrll]. 

We leave easy details to the reader. □ 



4.2 Estimates for ionization energy 

Theorem 4.2.1. 21 ) Let condition (3.5.2) be fulfilled and let N > Z - C Z^ . 

Then 

(i) In the framework of the fixed nuclei model 

(4.2.1) l* w <CZi. 

(ii) In the framework of the free nuclei model with N > Z — Q)Z^(Z~ <5 +aZ s ) 

(4.2.2) I* := E* N _, - E* N _, < zl (Z~ 5 ' + (aZ) 5 ') . 

Proof. Recall that Theorem 24.5.3 of [Ivrll] was proven simultaneously with 
Theorem 24.5.2; we follow the same scheme here picking up A = A^- Thus 
here and in the first part of the proof of Theorem 4.2.2 we estimate from 
above \* N {A). 

Again easy details a left to the reader. □ 

Theorem 4.2.2. 22 ) Let condition (3.5.2) be fulfilled and let N < Z - C Z 5 i . 
Then in the framework of fixed nuclei model under assumption (24.6.2) of 
[Ivrll] 



17 

18 



(4.2.3) {\* N + v)+<C{Z-N)*Z 



21 ) cf. Theorem 24.5.3 of [Ivrll]. 

22 ) cf. Theorem 24.6.3 of [Ivrll]. 




as a < Z 3 , 
as a > Z~3. 
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Proof. To estimate 1^ + v from above we follow exactly the arguments of 
Subsection 24.6.1 of [Ivrll] to estimate \n(A n ) > \* N . □ 

Problem 4.2.3. To prove the same estimate for (\* N + z/)_. 

Remark 4-2.4- To estimate \* N + v from below we pick up A = An-i] then 

\ N {A) = E N ^(A) - E N (A) < E* w _ x - E* w = \* N 

and we should follow the arguments of Subsection 24.6.2 of [Ivrll]. However 
in contrast to all other proofs of this Section here we should use spectral 
properties of Ha,v (or at least an estimate from above for its lowest eigenvalue 
after localization to supp# while in all other results we need an estimate 
from below for the same lowest eigenvalue after localization to supp6*) and 
to do so we need some uniform (i.e. with constants which do not depend on 
N) smoothness estimates for /4/v_i where /4/v_i is the minimizer for E/v_i(/4) 
as defined in Sections 3 and here (rather than as defined in Sections 1-2). 
While (5.A.1) is an analogue of (1.1.14), it is still not the same and while 
it implies some estimate it is not even remotely as good as we achieved in 
Sections 1 and 2. 

Sure p\\i is not very smooth either but it close to rather smooth p TF ; on 
the other hand minimizer A^ is an almost-minimizer for the one-particle 
trace problem studied Sections 1-2 but we don't know how close it to the 
minimizer (or one of the minimizers) of the latter problem. 

4.3 Free nuclei model: Excessive positive 
charge 

Theorem 4.3.1. 23 ) Let condition (3.5.2) be fulfilled. Then m the frame- 
work of free nuclei model with M > 2 the stable molecule does not exist 
unless 

(4.3.1) Z- N< Z b i(Z- s + (aZ) 5 ). 

Proof. Again we just repeat the proof of Theorem 24.6.4 of [Ivrll]. □ 
23 ) cf. Theorem 24.6.4 of [Ivrll]. 
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5 Appendices 



5. A Minimizers and ground states 

First establish a conditional existence of the minimizer 24 ) and the corre- 
sponding ground state of the original problem: 

Theorem 5.A.I. Let aZ < k* and let E* N < E* N _ V Then there exist a min- 
imizer A = An for the original multiparticle problem and the corresponding 
ground state U'/y. 

Proof. We know that as aZ < k* (with k* > which does not depend on Z 
or positions of the nuclei) E*(/4) is bounded from below; then ||V x A'\\ 2 is 
bounded from above for near-minimizer A' (but constants do depend on Z 
and (k* - olZ) here). On the other hand, if A' E and E N (A') < E N -i(A') 
there exists a ground state ^n(A'). 

Therefore if Ar^) £ ^o°° ^ s a minimizing sequence for E/y(/4) we have also a 
sequence W/v^/c)) with ||V/v(^(/o)|| = 1- Going if necessary to subsequence 
we can assume that A^) converges weakly in M' 1 and strongly in ££\ OQ for 
any p < 6; let A be its limit. 

One can prove easily that ^(^(k)) converge weakly in M' 1 and strongly 
in ££ 2 to \|/ and 

(H AV V, M/) = \\mJH A{k>v V N (A w ), V N (A W )) 

and then 

(H AV v|/, M/) + - = lim E A/ (^ ( , ) ) 

which is since A^) is a minimizing sequence and then W must be a ground 
state. □ 

Now in this framework we establish properties of these minimizer and 
the ground state: 

Proposition 5. A. 2. 25 - ) Let v|/ = N and A = An be a ground state and 
minimizer with energy E* N < E* N _ 1 - 

24 ) We do not know if it is unique. 

25 ) cf. Proposition 25.A.7 of [Ivrll]. 
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(i) W etf 1 andW = 0(e e l^l) as \x\ ->■ oo; 

(ii) Aetf 1 and A= 0(|x|" 2 ), V x A = 0(|x|- 3 ) as |x| -> oo; 

(Hi) Let N < Z. Then V* - V e tf 1 and = (Z - A/)^^ 1 + 0(|x| 2 ), 
V\4 = (Z - N)\x\- 2 + 0(|x|- 3 ) as |x| -»■ oo. 

Proof. Obvious proof using also an equation 
(5.A.1) -AA- = 

-2A/tr y V f (x,x 2 , ... ,x N )aj{D kx - A k {x))v k m{x, x 2 , ... ,x N )dx 2 - --dx N , 

is left to the reader. This equation is similar to (1.1.14) and is also derived 
from variational principles, its right-hand expression is 14 where A is the 
lowest eigenvalue of H^ y on Fock' space. 

5.B Zhislin' theorem 

Theorem 5.B.1 (Zhislin's theorem). 26) E^ +1 <E* N as N < Z. 

Proof. An easy proof repeating with obvious modifications proof of Theo- 
rem 25. A. 8 of [Ivrll] is left to the reader. □ 

5.C L. Erdos— J. P. Solovej' lemma 

We reproduce here Lemma 2.1 from L. Erdos, J. P. Solovej [ES3]. 

Lemma 5.C.I. There is a positive universal constant k* such that for any 
Z, a with Za < n* we have 

infinf H AV > - CZh 1/2 - zh~ 2 

N A 

2 

if CZ 3 < 5 < Ci with a sufficiently large constant C. 
26 > cf. Theorem 25.A.8 of [Ivrll]. 
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Proof. Consider pair of smooth functions 6q and 9\, 



6 2 + 9 2 = 1, supp^i C 6(0, 2r), 0i = 1 on 6(0, r), |V0 O |, |V0i | < CK 1 
with r = 5Z~^. 

Let xo be a smooth cutoff function supported on 8(0, 3r) such that 
|V%o| < Cr _1 and xo = 1 on 6(0, 2r). Let /4 be an average of A over 6(0, 3r). 
Let A) := {A - A)x , 6 := V x /A ; then V ® A = xoV <g> A + (A - A) <g> Vxo- 
Clearly 

/ B 2 dx< [ \V®A \ 2 dx<2[ x 2 o\V®A\ 2 dx+Cr- 2 [ {A-A) 2 dx 

</R 3 </R 3 </R 3 J B(0,Zr) 



< Ci / |V<gM 



dx 



for some universal constant Ci, where in the last step we used the Poincare 
inequality. Let y? be a real phase such that Vy? = A. Since xo = 1 on the 
support of 0i, we have 

0iH A , o 9i = e 1 e-'^H A _- A ^6 1 = e^H^ve^Or. 

After these localizations, we have 



(5.C.1) Hl z . A :-- 

N 



+ - 

J a </B(0,3r) 



^[^i(Wao -ft- (|V0 o | 2 + |V0i| 2 )0i 

j=l lx ' 

A/ 

>I]^ie"' V («A,.o-^M)e , '^i 



|V ®A\ 2 dx 
1 



2Cict 



e 2 dx 



with 



W(x) = 



+ Cr~ 



l(|x| < 2r) 



where 1(X) is a characteristic function of X. 
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Now we use the "running energy scale" argument in E. Lieb, M. Loss, M. 
and J. Solovej [LLS]. 



r i f°° 

(5.C.2) Y\\e 1 e-^[H A , tQ -W]e^e 1 >-/ N_ e (H A ,, - W) de 
. =1 L ij Jo 

>-£ N_ e (H^, -W)de- J N (^H^,o - W + e) de 

f°° e e 2 \ 

>-/ U_ e {H A , tQ -W)de- N (H A , i0 --W+-)de, 
Jo h V> V' 

where N_ e (/-/) denotes the number of eigenvalues of a self-adjoint operator 
H below — e. 

In the first term we use the bound H Ao0 > (D — A ) 2 — \B \ and the CLR 
(i.e. Cwikel-Lieb-Rozenblum) bound: 

(5.C.3) f 1 N- e {H Ao ,o - W) de < C T de [ (1/1/ + \B \ - e)f dx 
Jo Jo Jr 3 

<C I de [ {W - e/2)\dx+ C [ de [ {\B \ - e 2 /2^)f dx 
Jo Jr 3 Jo Jr 3 

<C Widx+Cfi* Sq dx = CZ^r^ + Cr~ 2 + CpS I B 2 dx. 



In the second term of (5.C.2) we use 

p 1 ?p7 
H Ao ,o ~ -W > - [(D - A f - — l(|x| < 2r)] 

+ ±{D-A ) 2 -\B \-^l{\x\<2r), 

and that 

( D -, ) 2 -^i( M < 2r )>( D -^-^>-(^ 

i.e. 

H Ao , ~-W>\(D- A f - 2(— ) 2 - |flb| - — 2 l(|x| < 2r). 
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Let (j, = 4Z 2 , then using Ce/fir 2 < e 2 /4fi for fj, < e (i.e. C < (5Z 2/3 ) 2 ), we 
get 

r°° e e 2 

(5.C.4) / M (H Ao , --W + -)de 

roc 1 J2 pfj, p 

< M (-(D-A ) 2 -\B \ + —)de<C de / (|6 | - e 2 /4^)f dx 

< / B 2 dx. 

Note that if Za < k* with some sufficiently small universal constant k,*, then 
(5.C.4) can be controlled by the corresponding term in (5.C.1). Combining 
the estimates (5.C.1), (5.C.2), (5.C.3) and (5.C.4) we obtain 

H AV > -czhi - cr 2 

and Lemma follows. □ 
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